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Abstract 

I Orthogonal Laurent polynomials in the unit circle and the theory of Toda-like integrable systems are 

■ connected using the Gauss-Borel factorization of a Cantero-Moral- Velazquez moment matrix, which is 
constructed in terms of a complex quasi-definite measure supported in the unit circle. The factorization 
of the moment matrix leads to orthogonal Laurent polynomials in the unit circle and the correspond- 
ing second kind functions. Jacobi operators, 5-term recursion relations and Christoffel-Darboux kernels, 
projecting to particular spaces of truncated Laurent polynomials, and corresponding Christoffel-Darboux 
formulae are obtained within this point of view in a completely algebraic way. Cantero-Moral- Velazquez 
sequence of Laurent monomials is generalized and recursion relations, Christoffel-Darboux kernels, pro- 

^ I jecting to general spaces of truncated Laurent polynomials and corresponding Christoffel-Darboux for- 

2^ ■ mulae are found in this extended context. Continuous deformations of the moment matrix are introduced 

and is shown how they induce a time dependant orthogonality problem related to a Toda-type integrable 

■ system, which is connected with the well known Toeplitz lattice. Using the classical integrability theory 
tools the Lax and Zakharov-Shabat equations are obtained. The dynamical system associated with the 
coefficients of the orthogonal Laurent polynomials is explicitly derived and compared with the classical 

I Toeplitz lattice dynamical system for the Verblunsky coefficients of Szego polynomials for a positive mea- 

sure. Discrete flows are introduced and related to Darboux transformations. Finally, the representation of 
the orthogonal Laurent polynomials (and its second kind functions), using the formalism of Miwa shifts, 
in terms of r-functions is presented and bilinear equations are derived. 
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1 Introduction 

This paper puts focus on orthogonal Laurent polynomials in the unit circle (OLPUC), a subject with strong 
links to that of orthogonal polynomials in the unit circle (OPUC), see [61], |56j and [57j. Despite that is well 
established that this matter is as source of interesting problems and applications in approximation theory, 
we are mainly interested in its connections with the general theory of integrable systems. 

Let us introduce here some notation that will be used along this article. We will denote the unit circle by 
T := {z G C : |z| = 1} , D := {z G C : l^l < 1} stands for the unit disk and ^[p,q] '■= spanjz"^, z~'P~^^, . . . , z'^} 
denotes the linear space of complex Laurent polynomials with their corresponding restrictions on their 
degrees while A[oo] for the infinite set of Laurent polynomials. When z € T we will use the parametrization 
z = e'^ with 6 € [0,27r). 

A complex Borel measure ^ supported in T is said to be positive definite if it maps measurable sets 
into non-negative numbers. When the measure ^ is absolutely continuous with respect to the Lebesgue 
measure of the circle d^, according to the Radon-Nikodym theorem, it can be always expressed using a 
complex weight function w, so that dfj,{9) = w{6)(i6. If in addition the measure is positive definite then the 
weight function should be a non-negative function on T. For notational simplicity we will use, whenever it 
is convenient, the complex notation d/x(z) = ie'^d/x(0). If ^ is a positive Borel measure supported in T, then 
the OPUC or Szego polynomials are monic polynomials P„ of degree equal or less than n that satisfy the 
following system of equations, the orthogonality relations, 

[ Pniz)z-''dfi{z)=0, k = 0,l,...,n-l. (1) 

It is well known the deep connections between orthogonal polynomials in the real line (OPRL) in [—1, 1] and 
OPUC (e.g. [3ll|T7j). Let us observe that for this analysis the use of spectral theory techniques requires 
the study of the operator of multiplication by z. The study of the matrix associated to this operator leads 
to recurrence laws. Both, OPRL and OPUC have recurrence laws, but the main difference is that in the 
real case three term recurrence laws provide a tridiagonal matrix, the so called Jacobi operator, while in 
the circle case one is lead to a Hessenberg matrix [37], leading to a more involved scenario to deal with 
than the Jacobi case (as it is not a sparse matrix with a finite number of non vanishing diagonals). More 
precisely, the study of the recurrence relations for the OPUC requires the definition of the reciprocal or 
reverse Szego polynomials Pi'{z) := z^Pi{z~^) and the reflection or Verblunsk}0 coefficients ai := Pi{0). 

^Schur parameters is another usual name. The definition is not unique and a; := — Pi+i(0) is another common definition. 
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With these elements the recursion relations for the Szego polynomials can be written as 



Pi 
p* 



) 



z ai 
zai 1 




) 



(2) 



There has been a relevant number of studies on the zeroes of the OPUC, see for example [10\ [T^ I15j . or 
[32l [35l SZl EI] , which have interesting applications to signal analysis theory, see [391 HD E2l [53] • Despite of 
that the situation is still far from the corresponding state of the art in the OPRL context. A second important 
issue is the fact that the set of Szego polynomials is in general not dense in the Hilbert space L^(T,/x). As 
it follows from Szego's theorem it holds that for a nontrivial probability measure on d/i with Verblunsky 
coefficients {a„}J^g the Szego's polynomials are dense in L'^{T, fi) if and only if ri^oC-*^ ~ l*^")!^) = 0- If d/x 
is an absolutely continuous probability measure then the Kolmogorov's density theorem follows: polynomials 
are dense in L^(T, /x) if and only if the so called Szego's condition Jj\og{w{9)d9 = — oo is fulfilled [60] . 

Orthogonal Laurent polynomials in the real line (OLPRL), where introduced in [421 143| in the context 
of the strong Stieltjes moment problem, i.e finding a positive measure /x such that its moments 



are known in advance. When the moment problem has a solution, there exist polynomials {Qn} such that 



which are called Laurent polynomials or L-polynomials. The theory of Laurent polynomials on the real line 
was developed in parallel with the theory of orthogonal polynomials, see \12\ [29l HQ] and [50]. Orthogonal 
Laurent polynomials theory was carried from the real line to the circle [62] and subsequent works broadened 
the matter (e.g. [251 |2T1 |23l |2l] ) , treating subjects like recursion relations, Favard's theorem, quadrature 
problems, and Christoffel-Darboux formulae. 

The analysis of OLFUC, and specially the use of the Cantero-Moral- Velazquez (CMV) |21j representa- 
tion is very helpful in the study of a number properties of Szego polynomials. Different reasons support this 
statement, for example as we mentioned while the OLFUC are always dense in L^(T, /x) this is not true in 
general for the OFUC, see |18] and [25], and also the bijection between OLFUC in the CMV representation 
and the ordinary Szego polynomials allows to replace the complicated recurrence relations with a five-term 
recurrence relation more alike to the structure of the OFRL. The representation of the operator of multi- 
plication by z is much more natural using CMV matrices than using Hessenberg matrices. This is the main 
motivation for us in order to take CMV matrices as a essential element in our scheme. Other papers have 
reviewed and broadened the study of CMV matrices, see for example [581 144] . Alternative or generic orders 
in the base used to span A[oo] can be found in [24] 

The approach to the integrable hierarchies that we use here is based on the Gauss-Borel factorization. 
The seminal paper of M. Sato [55] and further developments performed by the Kyoto school [26]- [28] settled 
the Lie-group theoretical description of the integrable hierarchies. It was M. Mulase, in the key paper |49) . 
the one who made the connection between factorization problems, dressing procedures and integrability. In 
this context, K. Ueno and T. Takasaki [63j performed an analysis of the Toda-type hierarchies and their 
soliton-like solutions. In a series of papers, M. Adler and F. van Moerbeke |3]-|9], made clear the connection 
between the Lie-group factorization, applied to Toda-type hierarchies (what they call discrete K-F) and the 
Gauss-Borel factorization applied to a moment matrix that comes from an orthogonality problems; thus, 
the corresponding orthogonal polynomials are closely related to specific solutions of the integrable hierarchy. 
See [45] and for further developments in relation with the factorization problem, multicomponent Toda 
lattices and generalized orthogonality. In the paper |8] a profound study of the OFUC and the Toda-type 
associated lattice, called the Toeplitz lattice (TL), was performed. A relevant reduction of the equations 
of the TL has been found by L. Golinskii p6] in the context of Schur flows when the measure is invariant 




(3) 




j = 0, . . . ,n - 1, 



(4) 
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under conjugation, another interesting paper on this subject is [48]. The Toephtz Lattice has been proved 
equivalent to the Ablowitz-Ladik lattice (ALL), [H [2], and that work has been generalized to the link 
between matrix orthogonal polynomials and the non-Abelian ALL in [20]. Both of them have to deal with 
the Hessenberg operator for the multiplication by z. 

Our aim is to explore the connection between Toda-type integrable systems and orthogonality in the 
circle from a different point of view. As we proof in this paper the CMV representation is a bridge between 
the factorization techniques used in [12] and the circular case. We will see that many results obtained in 
|12j about Christoffel-Darboux (CD) formulae, continuous and discrete deformations, and r-function theory 
can be extended to the circular case under the suitable choice of moment matrices and shift operators. 

Let us recall the reader that measures and linear functionals are closely connected; given a linear func- 
tional ^ on Ajoo] we define the corresponding moments of ^ as c„ := J£\z^\ for all the possible integer 
values of n € Z. The functional is said to be Hermitian whenever c_„ = c^, Vn € Z. Moreover, the 
functional is defined as quasi-definite (positive definite) when the principal submatrices of the Toeplitz 
moment matrix (Aj^), Ajj- := Cj_j, associated to the sequence is non-singular (positive definite), i.e. 
Vn G Z,A„ := det(cj_j)[^j-^Q 7^ 0(> 0). Some aspects on quasi-definite functionals and their perturbations 
are studied in [131 ttl] ■ 

It is known |34j that when the linear functional J£ is Hermitian and positive definite there exist a finite 
positive Borel measure with a support lying on T such that ^[/] = fj fdfi, V/ G ^[oo]- In addition, an 
Hermitian positive definite linear functional ^ defines a sesquilinear form ■ A^oo] ^ ^[00] 1— C as 

if J 9)^ = -^[/^li V/ifi^ £ ^[00] • Two Laurent polynomials {f,g} C A[oo] are said to be orthogonal with 
respect to .if if (/, g)^ = 0. From the properties of .if is easy to see that (•, is a scalar product and if fi 
is the positive finite Borel measure associated to .if we are lead to the corresponding Hilbert space L'^iT, /i), 
the closure of A^^y As we mentioned before {Pi}'^q denotes the set of monic orthogonal polynomials, 
degP; < for a positive measure fi satisfying ([1]) and therefore {Pi}1^q is an orthogonal basis of the space 
of truncated polynomials A[o^q] . 

In this work we allow for a more general setting assuming that .if is just quasi-definite, which is associated 
to a corresponding quasi-definite complex measure /i, see [33]. As before a sesquilinear form (•, •) ^ is defined 
for any such linear functional thus, we just have the linearity (in the first entry) and skew-linearity (in 
the second entry) properties. However, we have no symmetry allowing the interchange of the two arguments. 
We formally broaden the notion of orthogonality and say that / is orthogonal to g if {f,g)j^ = 0, but we 
must be careful as in this more general situation it could happen that {f,g)^ = but {g, /)^ 7^ 0. 

The layout of this paper goes as follows. In f|2]we present the application of the Gauss-Borel factorization 
of a CMV moment matrix to the construction of OLPUC, associated second kind functions, 5-term recursion 
relations and CD formulae. In ^ we perform a similar work using a more general sequence that the one 
used in |21j . This allows us to study snake-shaped (as are denoted in |24] ) recursion formulae. Moreover, 
the CD formulae we derive for these extended cases is the kernel for the orthogonal projection to the general 
space of truncated Laurent polynomials, A[p g] with p,q £ N. This is a large generalization of the CMV 
situation as the possible spaces of truncated Laurent polynomials is very particular, namely either Aj; or 
with / E N. To conclude, in §3] we study the deformations of the moment matrices to obtain the 
integrable equations associated, the representation of the OLPUC and its associated second kind functions 
using r-functions and corresponding bilinear identities. 

2 Orthogonal Laurent polynomials in the circle, LU factorization and 
the CMV ordering 

In this section we use the Gauss-Borel, also known as LU or Gaussian, factorization of an infinite dimensional 
matrix, that we refer to as moment matrix, to derive bi-orthogonal Laurent polynomials in the unit circle 
(BOLPUC) to the given measure /j, -that, as we will see, are closely related to the Szego polynomials-, 
associated second kind functions in terms of the Fourier series of the measure, their recursion relations and 
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corresponding CD formulae very much in the spirit of [12]. The key idea is to use the results of |i21j| to 
construct a very specific moment matrix that will lead to the mentioned results. 

2.1 Biorthogonal Laurent polynomials 

We are ready to, using a CMV moment type matrix, find a set of bi-orthogonal Laurent polynomials, 
its connection with Szego polynomials and corresponding determinantal expressions. For this aim, let us 
consider the basic object fixing the CMV order of the Fourier family {z-'ljgz. This order allows us to work 
in realm of the semi- infinite matrices avoiding in this way the less convenient scenario of bi- infinite matrices, 
very much as in the multiple OPRL situation [12j . 

Definition 1. We denote 

Xi(^) := (l,0,z,0,z2,0,...)T, ;t2(^) := (0,1,0, z, 0,^2,... )T, xliz) ■= z-^Xaiz~^), a = 1,2. 



Xiz) ■■= Xi{z) + xliz) = ,z,z 



-2 



X*iz) := xliz) + X2iz) = iz'\ 1, z~\ z,...) 



With these sequences at hand and with a given quasi-definite complex Borel measure /U supported on T 
we define the CMV moment matrix 



Definition 2. The CMV moment matrix is the following semi-infinite complex-valued matri^ 



9 ■-- 



X{z)x{z)'^ ^Im{z) 



(5) 



The LU factorization, which will play an important role in what follows, is 

g = S^^S2, (6) 
where 5*1 is a normalizecH lower triangular matrix and 5*2 is an upper triangular matrix. Hence we write 



(I ...\ 

(5i)20 ('S'l)21 1 



v 



So 



( (5'2)oo (5'2)oi (5'2)o2 • • 
(^2)11 (^2)12 .. 
(^2)22 



V 



This Gaussian factorization for the moment matrix makes sense if all the principal minors are non-singular, 
which is precisely the quasi-definiteness condition for the measure [i. More on the algebraic Gauss-Borel 
(or LC/) factorization and it's connection with integrable systems can be read in j30] . 
With the aid of these matrices we consider the sequences 



^1,1 := Sxx\, ^1,2 := 5'iX2, 

which can be written as semi-infinite vectors 



$2,1 := (52"')^Xi 



^2,2 := (52"')^X;, 



^lAz) 



{^ti.z)\ 
^S(^) 



^2,a(^) 



4!i(^) 



V 



^The reader should notice that if is a positive measure then g is a definite positive Hermitian matrix; i.e., g = ■ 
^The coefficients in the main diagonal are equal to the unity. 
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for a = 1,2. The corresponding components (p^\ are polynomials of degree / in variable z, while (^['2 
polynomials of degree I + 1 in the variable which vanish at 2; = 00. Inspired by the multiple orthogonal 
case |12] we define the following sequences of Laurent polynomials 

$i(z) : = Six{z) = $1,1 + «>i,2, ^z) : = {S^yxiz) = ^2,1 + ^2,2, 

which are semi-infinite vectors that we write in the form 



^l(^) 



\ 



J 



V 



with its coefficients being Laurent polynomials. 

As we said in blithe measure fi has an associated sesquilinear form (•, ■)_^ acting on any pair of Laurent 
polynomials in T, ipi{z) and (p2{z), as 



if 1,^2): 



ipi{z)ip2{z ^)dfi{z). 



(7) 



From the definition is clear that g (whose principal minors should not vanish) is the matrix associated 
to The reader can check that the principal minors of g are exactly the Toeplitz minors A„ (one 

matrix is obtained from the other using permutations). As we are going to use quasi-definite measures the 
factorization condition will hold in the subsequent work. 

We recall the reader that given two linear spaces V, V and a sesquilinear form 

(•,.)^: VxV ^ C 

{x,y) ^ {x,y)cf 

we say that sets X CV and Y cV are bi-orthogonal if {x,y)y = for all x G X and y gY. 

Theorem 1. The sets of Laurent polynomials {ipi^}'^^ and {(/?2^}^q are bi-orthogonal in the unit circle 
with respect to the sesquilinear form defined in ([7]), that is 



I (0 (fc)\ 



^f\z)^^^'{z-')df,{z) = 5i„ 



-Jk) 



l,k = 0,1,... 



(8) 



Proof. We compute 



^l{z)<i>2iz-^)'dfi{z 



<i>i{z)'^2{z)^dfi{z) 
X(.z)x{z)^ dfi{z) 



Si 
I. 



Orthogonality relations ^ can alternatively be expressed as follows 



"^2 



□ 



{^r,z')^ 
'r'\z')j. 



(pf)(z)z-«d^(z) = 0. 



/ fc (2^+1) \ 



k — I ^ . . . ^ I 1 5 

V9f+^^(Z)Z-M^(z) = 0, k-- 

(^f^(z-i)z'=d/u(z) = 0, k-- 

T 

ipf^^^\z-^)z''dfi{z) =0, k = -l,. 



I, . . . ,1, 
I,. ..,1-1, 
J. 



(9) 
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Proposition 1. Given a positive definite measure fi there exists /i; € M, / = 0, 1, . . . , such that 

ipy = hi ip\' 

Proof. See Appendix [XI □ 

Therefore, are proportional Laurent polynomials and bi-orthogonality ([5]) implies that {Vi^l^^o ^^'^ 
are sets of orthogonal Laurent polynomials for the positive measure fi, that is 

{^?,vf^)^= <f ip?iz)4'\z-^)df,iz) = 5i,khi l,k = 0,l,... 



(10) 

{ip^\ip^i^)^= f ip^j\z)4'\z-')dfi{z) = 6i,khr' l,k = 0,1,... 

We are now ready to write the relationship between these CMV Laurent polynomials and the Szego 
polynomials Pi introduced previously. 

Proposition 2. // the measure fi is positive definite we have the following identifications between the CMV 
Laurent polynomials, the Szego polynomials and their reciprocals 

^f\z) = Z-'P21{Z), ip?^'\z) = Z-'-lP2Vi(z). (11) 

Proof. To check it just observe that 

zVf iz)z-''dn{z) = 0, k = 0,...,2l-l. 

T 

Hence, z^ip^i^\z) has the same orthogonality relations that P21 and both are monic polynomials of degree 
21; uniqueness leads to their identification. In a similar way we proceed for the odd polynomials. Indeed, 

^i+i^(2/+i) (^)^-fe^^(^) ^0, A; = 1, . . . , 2/ + 1, 

T 

that is, z^^^ ip^''^^\z) has the same orthogonality relations that the polynomial -P2i+i (th^it makes them 
proportional) and both are equal to 1 at z = 0, consequently they are the same. □ 

Using the Verblunsky coefficient^ we can write 



^f'\z) = a2iz-' + --- + z', 
^f'^'\z) = z~^-' + ... + a,i^,zK 

For later use, and in addition to the reflection coefficients, it is also useful to define the sequence pi :- 
y^l — related to {hi}'^^ by 

Pi 



2 _ 



hi-i 

valid for I > 0, with po = 0. 

In the general quasi-definite case, we can perform a very similar construction. We write 



(20, 




= o^z-^ + - 




(2/+1), 






^ "2/4 


(20, 
f2 1 






+ ••• H 


(2/+1), 
^2 1 




- h-^ 

~ "'21+1^ 


+ . . . ^ 



(12) 

iz\ 

f'^^ z^ ^''^ 

"'2/+l"2Z+l^ ' 



*See the Introduction. 
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and also pi := 0, pf := t^J = 1, 2, 



(2) ± 

= ce^^''. The reason for the notation stands in the following fact. Given a quasi-definite measure fi we 
can find two families of monic orthogonal polynomials such that 

Pl^\z)z-''dp{z) = A; = 0, 1, 1, 

(14) 

z''pP{z-^)dp{z) = fe = 0,l,...,/-l. 

JT 

If we call = P^^^^(O), and = pI^'^\o) then these coefficients are related with the Laurent Polynomial 
coefficients and we can obtain the quasi-definite version of Proposition [21 that is 



„(i) 



hi 



The Hermitian case can be considered a particular reduction where 



^r\z) = z-^ps\z), ^r'\z) = z-^^^p^iz), 

^f\z) = -h,lz-'pS\z), ^f^'\z) = h,l^,z-'-^P!il\{z). 



{21+1), 



(15) 



In what follows := Yli~'j=o9iJ-^iJ denotes the / x / truncated moment matrix and x^^^ is the truncated 
vector consisting of the I first components of x- With this notation we can express the bi-orthogonal Laurent 
polynomials in different ways: 



Proposition 3. The following expressions hold true 



{S2)ii{0 



det oH 



det 



50,0 50,1 
51,0 51,1 



1) (5"+^l)-V['+^] 



90,1-1 
91,1-1 



\ 9l,0 91,1 ■ ■ ■ 9l,l-i 



and 



4'^(--) = (52),H(x«)^-(/l)n5W)-^ 



/ 9o,i \ 
9i,l 

\9l-i,lJ 







(^[m])t(^[m])-i 





VV 

/ 50,0 
51,0 



det ^I'+i] 



det 



50,1 
51,1 



5/~i,o 5«-i,i 



5o,« \ 
9i,i 



9i-i,i 



WHz) x«(^) ... x«(^)/ 



Proof. See Appendix lAl 



I > 1, 



(16) 
(17) 

(18) 



(19) 



(20) 



I > 1. 



(21) 



□ 



Similar expressions hold for replacing x by xi; a = 1; 2, and for ip^^^ replacing x by = !> 2 



(0 
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2.2 Second kind functions 



In this subsection we introduce the second kind functions associated with the orthogonal Laurent polynomials 
discussed before. First we present determinantal expressions, then we connect them with the Fourier series 
of the measure and also with corresponding Cauchy/Caratheorody transforms. 

Definition 3. The partial second kind sequences are given by 

Ci,i{z) := {S^^Yxliz), Ci^z) := {S^')^X2iz), C2,i(z) := S^xli^), C2,2{z) := S2X2{z). 
and the second kind sequences 

Ci{z) := (Sri)tx*(z), C2{z) := S2X*{z). 

Observe that 

Cl = Cl,l + Ci^2i C2 = C2,l + C2,2- (22) 

We have the expressions as semi-infinite vectors, being its coefficients what we call second kind functions 



Ci,a{z) :-- 



\ 



lcfl{z)\ 

ciiiz) 



J 



V 



Ci{z) :-- 



fc['Hz)\ 

C['\z) 



J 



V 



C2iz) :-- 



I 



(Cf\z)\ 

C^\z) 

V : / 



The expressions just given for the second kind functions are, in principle, formal, as the matrix products 
lead to series, not necessarily convergent, instead of finite sums. In fact (as we will show in Proposition 
[5]) they are well defined in terms of the bi-orthogonal Laurent polynomials and truncated Fourier series of 
the measure, with convergence ensured in some annulus centered at the origin of the complex plane. The 
coefficients of these sequences are called second kind functions. 

One can find analogous determinantal type expressions as those for the OLPUC given in Proposition [3l 
if we define 



r(0 .. 



(fc) 



k>l k>l 

Proposition 4. The second kind functions have the following determinantal expressions for I > 1 

( 90,0 50,1 ••• 9o,i \ 
51,0 51,1 ••• 9i,l 



C?iz) 



det ^I'+i] 



det 



5/-i,o 9i-i,i 



9i-i,i 



\^r?i(z) Tl{{z) . . . Tl'iz) J 

\ 



.(0 



(23) 



and 



cS\z) 



1 



det 5(W 



det 



50,0 50,1 
51,0 51,1 



\ 91,0 91,1 



9o,i-i 
9i,i-i 

9i,i-i 



rS(--) 

rS:](^) ; 



(24) 



Proof. See Appendix [Al 



□ 
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Now, we introduce the following definition 
Definition 4. 1. For the bi-orthogonal Laurent polynomials Lpf^ and ip2^ we use the notation ip^ {e^^ 

2. Let 

f-2TV 

2tt 



Cn = ^ e-'"^d/x(^), F^(n)= V c„«", 
Jo 



be the Fourier coefficients and the Fourier series of the measure fj,. 
3. For each integer k we introduce the following truncated Fourier series 

n>—k n<—k 

Observations 



1. It holds that Cn(/^) = C-n{^j), hence c_n = c„ for real measures. Consequently, 

F^i^ = F^^i^) = F,iz) = F,{z-'). 



2. The Fourier series always converges in D'(T), the space of distributions in the circle, so that J^^ F^{9)f{d)dd 
J^"" f{e)dij{B), V/ G V{T), here V{T) denotes the linear space of test functions on the circle. For an 
absolutely continuous measure dfi{9) = w{6)d9 we can write dfi{6) = F^{9)d9. 



3. We will also consider the Laurent series Ffj_{z) = J2'^=-oo ^nz"', for z G C. Notice that F^^^^+F^^ = F^. 

4. Let D{0;r,R) = {z € C : r < \z\ < R} denote the annulus around z = with interior and exterior radii 
r and R and R± := (limsup^^o^ \/|'^±"l)^^- Then, according to the Cauchy-Hadamard theorem, we 
have 

• The series F^^^(z) converges uniformly in any compact set K, K C D{0; 0, R+). 

• The series k^-^) converges uniformly in any compact set K, K C D{0; oo). 

• The series F^{z) converges uniformly in any compact set K, K C -D(0; 

5. According to F. and M. Riesz theorem if c„ = 0, n < 0, then /i is absolutely continuous with respect to 
the Lebesgue measure on T [54]. In fact if f{z) = ^ j^^ T^T^^' then w{6) = limr_j.i f{re^^) S L^{T) 

and d/i = w{9)d9; therefore, in this case we have f{z) = Jj ^^du, a holomorphic function in 
D. Moreover, the set of w{9) = /(e'^) G L^(T) with = 0, n < 0, is isometric to the set of 
holomorphic functions in B with limits when r ^ 1 in L^(T), observe that w = Yl'n=o Cn e"'^ = F^. 

Proposition 5. The partial second kind functions can be expressed as 

C[l = 27r ^ ^gz-'-'F^-2,_,{z), ii_<|.|<oo, 0^1 = 2^ ^^iy''F^^},_,{z), < \z\ < R^, 
Cg = 2vr 5; i?;^<N<oc, Ci% = 2n ^flz-'-'F^Jiz-'), < \z\ < R-_\ 

|fc|<oo |A:|-Coo 

and the second kind functions as 
Cf^ = 2T:^f{z^^)z-^F^{z), R^<\z\<R+, C^^'^ = 2Tr^f\z^^)z-^ F^{z-^), R^^ < \z\ < Rz\ (25) 



Here is used just to indicate that the sum is finite. 
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Proof. Prom the formal definition of Ca^j a,b = 1, 2, and the aid of the Gaussian factorization of the moment 
matrix g, we have 

Ci,i = (32^ i x{u)x{u)^Mu)x\{z), Ci,2 = (52-^)t / x{u)x{u)'^d{I(^)X2{z) 

C2,i = Sii x{u)x{'<j)'^^l^{u)Xi{z), C2,2 = Si f x{u)xiu)^(iKu)X2{z), 



We recall that ((S'2 ^)^x(^))^'^ = ¥'2^(^) ('S'ix(^t))^'^ = fi\u) and expand the matrix products involved, 
without any interchange of integrals and summation symbols, to get 



|fc|<oo 



n=0 

00 „27r 



|A.i<oo 



n=0 

00 „27r 



|fc|<oo n=0 "'O IfcKoo n=0 •'^ 



IfcKoo 



Using the Fourier coefficients c„ we write 

00 00 



|A;|<Coo 



n=0 

CO 



|fc|<oo 



n=0 

00 



eg = 2. 9^a(E«"""0' ^S = 2- E rfi(E-^— 

fc|<oo n=0 |fcl<oo 



n=0 



from where the desired result follows. 



Example. Assume that 



so that 



d/x(0) = w{e)de, 



w 



cosnO „„fl , . 
— = e + e'^"^ ^ cos (sm 6*) , 



n=0 



n! 



where the weight w is a smooth positive 27r-periodic function. Then 



2V E:=oV + e^ j fc>0, 



with annulus of convergence D = {z E C : < \z\ <1}. Notice that in this example 

F^{z)=e + 



e^ + e^ 



and therefore 



ci'\z) = 27rz-'^^\z-')F,{z), 



ci\z) = 27rz-'^'l>{z-')F^{z). 



(26) 



□ 
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Proposition 6. The formal series for ^(z) (where a = 1,2) defined in Proposition^^ can be expressed in 
terms of the Fourier series of fi and consequently are convergent in corresponding annulus in the complex 
plane. More precisely 

rg(z) = 2vrz-^(^-)-i(FW_,^(z-i) + R:,' < \z\ < Rz\ 



wher^ J{j) = being [p] the integer part of p. 

Proof. See Appendix El □ 

Notice that for / = we have 

rf^{z) = 27rz-^(j>iF^(z-i), R-' < \z\ < Rz\ 

rg(z) = 27rz-^(j')-iF^(z), R^ < \z\ < R+. 

Now we will justify the name we have given to these functions and show a Cauchy integral representation 
of them. We will prove it in two different scenarios depending on the measure. First, with a more measure 
theory taste, assuming that fi is positive and using the Lebesgue dominated convergence theorem. Second, 
with a more complex analysis taste, considering absolutely continuous complex measures d/i = w{6)d6 when 
tt; is a continuous complex function. 

Theorem 2. Assume a positive measure dfi{9) or a complex measure dii{9) = w{9)d9 with w a continuous 
function. Then, the second kind functions can be written as the following Cauchy integrals 

(0 = ,-1 / !f^dM^, c« = Z-' I ^^^d^.{u), \z\ > 1, 

Jj u- z ^ ' JfU- z ^ 

= -Z-' I ^^^dM^, = -Z-' I ^^^dMn), \z\ < 1. 

Jj u- Z ^ ' Jj u- z ^ 

Proof. From the definition of Ca^b ^-nd the aid of the Gaussian factorization of the moment matrix g, 

(5f i)t = {S^^)^g^ or 52 = S^g we get 

= (52-^)t j> x{u)x{uydi^)xl{z) =Y.[f ^2{u)u-^z-^-'Mu)), 

n OO „ 

Ci,2 = {S:,^ f x{u)x{nyMu)X2{z) = E ( / Mu)u^^'z"Mu)), 

/. oo „ 

C2,i = Sij> x{u)x{n)^Mu)x{{z) = E ( / ^l{u)u-^z-^-^d^i{u)), 

^2,2 = S,f x{u)xiu)^d^{u)x2{z) = f](f f i(n)n-+izM^(n)). 

For the series in these expressions we have: 

1. The series X^^o ^""-^""'^ converges uniformly in the u variable in any compact set iiT C {tt S C : 
\u\ > |^;|~"^} to {z — u^^)^^ and if \z\ > 1 then we can take K such that T C K. 

2. The series YlnLo u^'^^z^ converges uniformly in the u variable in any compact set K C {u ^ C : \u\ < 
\z\~^} to —{z — u~^)~^ and if \z\ < 1 then we can take K such that T C K. 



®The reader can check that z'''^' — X^''\^) 
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Let us assume a positive measure. The corresponding m-th partial sums are 



E 



1 - (nz)-i 



U- 



1 - (uz) 



n=0 ^ ' n=0 

If we write u = e'^ and z = \z\ e'^''^^ we have 

l-iuz)-""-^ 2 l-2[z|-('"+i)cos((m + l)(6l + argz)) + 



1 - (txz)-! 



1 — 2\z\ 1 cos(6' + arg z) + \z\ 



For [z| ""^ < 1 we have the fohowing inequahties 

0< l-2|z|-('"+i) cos((m + l)(0 + argz)) + Izl-^C^+i) < (i + |2[-(r«+i))2 < 4^ 
< 1 - 2|zl"^ cos(6l + argz) + \z\'^) > (1 - 

so that, for n G T, we infer 



1 - fuz)-i 



'"-I 2\z 

< 



-1 



1- 



Ul > 1. 



Similarly, we conclude that 



1 - 



< 



1 - \z\ 



1 — uz 

Thus, for n G T and a = 1,2, we have the control bounds 

m. „ 

\^^^\u)u--z—^ < 

n=0 



\z\ < 1. 



1 - z 



n=0 



< 



1 - \z\ 



\z\ > 1, 



\z\ < 1. 



Consequently, as the Laurent polynomials ipa are measurable functions in T, the Lebesgue dominated con- 
vergence theorem leads to the sated result. 

Finally, if we assume that d/i = w{9)d9, with w a continuous complex function; we can always write 
w{6)d6 = F{u)y^, u = e'^, with F a continuous function on T. Then, recalling the uniform convergence of 
the geometrical series involved and the fact that the Laurent polynomials ipa are continuous functions on T, 
we can interchange integral and series symbols arriving to the expressions 



1- lAEZo^aiu] 



u z 



n „— n— 1 



Mu) = i:n=o(k^a{u)u~-Z~--^' 



)d^{u) ] for \z\ > 1. 



2- /tt (En=oMu)u^-''z'')dKu) = En=0 ( fj ^a{u)u-+^ Z-)d^,{u))^Or \z\ < 1. 



□ 

The result motivates the name given to these functions [46]. These expressions can be also written as 
Geronimus transforms, for that aim we just need to recall that for u € T 



z — u 



2z\ 



-1 



u — z 
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and therefore for / > 1 

C7i;i = 1 / ^^±^^y.^\u)Mu), Cil = ^i±^v^S')(n)d^(^), \z\ > 1, 

ZZ U — Z ZZ J'Y U — Z 

Cfl = / ^i±^4')HdM^, eg = / ^i±^,,f)(n)d/u(n), \z\ < 1. 

U — Z ZZ 11 — Z 

For / = 0, we just obtain (up to constants) 

^(H+c(z-i)), 

where C{z) is the Caratheodory transform of the measure: 

1-217 i9 _|_ 

c{z) := / -^Me). 

Jo 6 —z 

Example. Let us assume that for \z\ > 1 we have 

and that the only singularity of F inside B can lay exclusively at z = 0. On the one hand, as F has a Laurent 
expansion of the form F{z) = X^feez -^fc-^*^ ^^"^ there no singularities different from z = we conclude that 
the coefficients of the Laurent series are just the Fourier coefficients of the measur^ 

On the other hand, as the Laurent polynomial ipi has its singularities at z = 0, and z~^ € D, then the set 
of singularities of the integrand that lay in D in (j27p is {0, z~^}. Hence, the residue theorem gives 

= 2..- (Res...-. (^fVm^) + Res„.o ( ^^M^ 
V \ u — z^^ / \ u — z'^ 

We now evaluate the residues, the first one may be computed noticing that ipf'' {u)F{u) is analytic at u = z^'^ 



._-.(5:iJ^)=rf>(.-'m.-') = c!' 



u — z 

where we have used Proposition [5j For the residue at n = we use the Laurent series around this point 

oo 



k j n=0 

SO that 



2vr.-Res..o (^^|^) = " E ^'IW""^' = " E E --n-i^-^ = -G 



k+j+n=—l k n=0 

n>0 



and we finally get 



C2,l — C2 — 6*2,2, 



in agreement with (j22p . 

The application of the residue theorem to the formulae in Theorem [2] leads to expressions for the second 
kind functions in terms of residues: 



^This property does not hold if F has different singularities from z = inside D, as we can not choose T as the integration 
circuit to apply the Cauchy formula for the coefficients of the Laurent series. 
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Proposition 7. Let us assume dfi{9) = F{u)f^ then: 



1. WhenF is an analytic function in C\ID) but for a set of isolated singularities, and if denote by {^j,— }j_q 
the set of distinct points obtained from the union o/zo,_ = and the set of singularities of F C\D, 
then, for z ^T, 



Ciliz) = 27TZ-^ ^'i\z~')F{z)d{\z\ - 1) - Res.=o (v?^'^(n)F(n-i))z + ^ ( ^^^^^ Res„^,-i F{u-' 



(0, 



(0, 



2. If F is an analytic function in D but for a set of isolated singularities, and denote by {^j,+}jio 

of different points obtained from the union of Zq = and the set of singularities of F at B. Then, for 
z^T, 

P+ if) ( \ 

Cg(z) = 27rz-i UP{z-^)F{z-^)e{\z\ - 1) - Res„=o {^^? {u)F{u))z + J] (^^^^ Res.=., F{u] 



z^+ — z 



ci%iz) = 27TZ-' 



ipf\z-')Fiz-yil - + Res„=o {^?iu)Fiu))z - ^ (^Li£Zi±^ Res,^^^ F{u] 



To conclude this section we give some summation rules that are derived using the geometrical series 
Proposition 8. The OLPUC and its corresponding partial second kind functions satisfy 

\z > \z\. 



/=0 



for a = 1,2. 

Proof. See Appendix lAl 



CO 








oo 








/=o 




oo 




EcS 




/=0 





z — z' 

1 



z — z' 



\z\ < \Z 



□ 



2.3 Recursion relations 

We are about to derive, using the Gaussian factorization, the CMV recursion relations and obtain in this 
way the well known CMV five diagonal Jacobi type matrix for the recursion of the Szego polynomials. Let 
us begin with the following 

Definition 5. Given the canonical basis for semi-infinite matrices Eij,i,j £ Z+, we define the projections 



III :- ^E2j,2j, 

j=0 



U2 :- ^£^2j+i,2j+i, 
i=o 
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and the matrices 



Ai 



,7=0 
oo 

3=0 



A2 •.= y^^Ei+2j,3+2j, 
3=0 

T:= Ai+A^ + Si,iA^ 



The matrix T, 



T 



/O 
_1_0 
0^ 

o_i 


0_0 
0^ 
0_0 
0^ 




V 



1 0]0 

g oio_g 
oji 
goio_g 
oTo 

OilO_0 

ojo 

OlO 1 



ojo 

OlO 
I 



0]0 

o_oig_g 

ojo 
g_oig_g 

1 oTo 

0_0l0 

oji 
o_oig_o 
o"|o 
no 



I 



is a central object in this paper, as its dressing -its orbit by conjugations- gives the pentadiagonal CMV 
Jacobi type matrix. 

It is immediate to check that 



Proposition 9. The following relations hold 

^iX{z) = zUixiz), Mx{z) = z-^U2x{z), 

^Jx{z) = {z-^^i - Eo,oA)x{z), Ajxiz) = {zU2 - Ei,,A^)x{z), 

'^Xiz) = zxiz), 
T^x{z) = z-\{z). 



(28) 

(29) 



With the aid of these conditions we characterize the moment matrix g as verifying a symmetry constraint, 
which we call string equation, from where the recursion as well as the CD formulae will be derived. The 
symmetry is detailed in the following 



Proposition 10. The CMV moment matrix fulfils the following condition 

Tg = gT. 

Proof. To prove the previous result we proceed as follows 

Tg = <f{Ai+Aj+Ei,iA^)x{z)x{z)''dfi{z) = I zx{z)x{z)'^ Mz) 
= i X{z){z~^x{z))^<^l^{z) 



(30) 



/ 

JT 

I x{z){{Al +A2 + E^,^A)x{z))^Mz) 

JT 



= g{Ai + A^ + A'^Eofi) 
= g{A, +Aj + Ei,iA^) 
= 5T. 



□ 
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We now proceed to dress the T matrix in two ways 
Definition 6. We use the notation 

Ji '.= SiTS^ ^, J2 S2TS2 ^■ 

The fohowing proposition is triviahy derived from the above definition and pO 
Proposition 11. 

Ji = J2. 



Consequently, we introduce the CMV Jacobi type matrix 

Definition 7. We define J := Ji = J2. 

The matrix J has a five diagonal structure; as easily follows when one observes that Ji has zero coefficients 
over the third upper-diagonal and that J2 has all its coefficients equal to zero under the third lower-diagonal. 
More specifically, the structure is 



J 



/ * * ] 1 

-|- * I * 








+ 1* 
~d" o|o 
_o_o 10 

OiO 



]o ]o ]o 

OlO OlO OlO 

*|i olo 

* I* 10 
-H- h.- 

* I * * I i 



oTo o' 



_0 

0" 





ojo 

10 



V 



I 



-I- 1 * * I * 

-4 \— 

OlO * I* 

ojo ±ll 

OlO OlO 

o]o o]o 



10 

o|o o' 

jO 

* il o' 



* I* * 
+ j* * 



/ 



where, * is a possibly non-vanishing term and -|- is a positive term. In fact, using the LU factorization 
problem we are able to completely characterize J in terms of the Verblunsky coefficients. 

Proposition 12. 1. The non-vanishing coefficients of J are 



J2k,2k-1 
J2k+l,2k-l 



2 (1) 

2 2 
P2fc+lP2fc) 



J2k,2k 
J2k+l,2k 



r,(2)„,(l) 



2 -(2) 
P2k+l'^2k ' 



'^2fc,2fc+l 
J2k+l,2k+l 



-a 



(1) 

2fc+2' 



J' 



2k,2k+2 



1, 



"«2fc+l«2fc+2' «^2fc+l,2fe+2 - a^k+r 



2. We have the recursion relations 



(2k) (2k+2) (1) (2k+l) 

M =n - «2fe+2'/'l 



-(2) (1) (2fc) 2 (1) (2fc-l) 
'^2kO'2k+m -P2k0^2k+m 



{2k+l) 
ZLp\ 



_(2) (2fc+2) _(2) (1) (2A:+1) , 2 -(2) (2fc) , 2 2 (2fc-l) 
"k+lV'l - a2fc+l«2yt+2V'l + /'2fc+ia2fcVl + /'2fe+l/'2feV'l 



(31) 
(32) 



Proof. 1. See Appendix 1X1 
2. We have 



(2fc) J (2fc+2) , J (2A:+1) , j (2k) , ^ (2fc-l) 

ZLp^ — J2k,2k+2Vl + J2k,2k+l^l + 'J2k,2kVl + -'2fc,2fc-lV'l > 

(2A;+1) J (2fc+2) , J (2fc+l) , j (2k) , ^ (2k-l) 

Zfi = J2k+l,2k+2'^i + J2k+l,2k+Wi + 'J2k+l,2k<fi + J2k+l,2k-l'Pi 



□ 
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As Ji$i = SiTS^^Sixiz) = z<^i, the sequences <^i,'J?2 have a five term recurrence formula. However, 
although there are five non- vanishing diagonals the recurrence relations do not have more than four non-zero 
terms, explicitly 



(2k) J {2fe+2) , J (2k+l) , J {2k) , J {2k-l) 

-= J2k,2k+2^i + J2k,2k+Wi + ■hk,2k'^\ + J2k,2k^Wi 



(2k+2) _ (1) (2fc+l) _ -(2) (1) {2k) _ 2 (1) {2k-l) 
^1 ^2k+2^1 "2fc "2fe+l'^l P2fc"2fc+l¥'l ' 



{2fc+l) J (2fc+2) , J (2fc+l) , J {2k) , J {2k-l) 

Zfi = J2k+l,2k+2fi + J2k+l,2k+lfi + J2k+l,2k^i + J2k+l,2k-Wi 

-(2) (2fc+2) _(2) (1) (2fc+l) , 2 -(2) {2k) , 2 2 {2k~l) 
= ^2k+m - "2fc+l«2fc+2'^l + te+ia^Vl + P2k+lP2kn 

to them we can add the truncated relations for A; = 0, 1 

(0) (2) (1) (1) (1) (0) 

(1) _(2) (2) _(2) (1) (1) , 2 (0) 

It is also possible to build recursion relations multiplying by , thus we have 
Proposition 13. The OLPUC have the following recursion relations 



(33) 



z 



-lA2k) _ (1) (2fc+l) ^(1) -(2) {2k) 2 (1) {2k~l) 2 2 ,,(2^-2) . .. 

— "2fc ~ ^2k "2fc+lV'l + P2k'^2k-lVl + P2k-lP2kV^l i l-^^j 

1 (2fc+l) _ (2fc+3) _ (2) (2fc+2) _ (1) (2) (2fc+l) _ 2 -(2) {2k) , . 
Z - fl '^2k+3'^l ^2k+1^2k+2^l P2k+l<^2k+2'Pl ' 

where we have to add the truncated relation 

-1 (0) (1) -(2) (0) 

Proof. Using that S'lT^S'j^-'^^i = SiT^ S^^ Sixiz) = z~^^i we need to calculate the coefficients of SiT^ S^^ = 
82'^'^ as we did with Ji, to obtain the desired result. □ 

With the previous result we can get 
Proposition 14. The coefficients pf verify the following relations 

pl = l-4^a^^\ k>0. 
Proof. See Appendix O □ 

The following results were found previously in [2Tj using an alternative derivation. Recalling the explicit 
form of the CMV Jacobi type matrix J provided in Proposition [12] we get 

Proposition 15. When the measure p is positive the recursion relations for the OLPUC can be expressed 
in terms of the Verblunsky coefficients as follows 



{2k) {2k+2) {2k+l) _ {2k) 2 (2fc-l) 

Zip\ ' = ip\ - a2k+2^\ - Oi2kOC2k+Wi - P2k'^2k+l^i 

{2k+l) _ (2fc+2) _ {2k+l) , 2 - (2A:) , 2 2 (2fc-l) 

M =Ol2k+l^\ - a2k+ia2k+2f l +p2k+l»2kVl +p2k+lP2kf l 

-1 {2k) {2k+l) _ {2k) , 2 {2k~l) , 2 2 (2A,-2) 

Z V^i =a2kn '-a2ka2k+W\ ' + p2ka2k^W\ + P2k~lP2k'Pl 

-1 {2k+l) (2A:+3) _ (2A:+2) _ {2k+l) 2 - {2k) 

Z V\ -a2fe+3</'i - a2fc+ia2fc+2V5l - P2k+l(^2k+2V\ , 



(36) 



(37) 



where k > 0. 

Using the Szego polynomials and their reciprocals recursion relations (|31|) and (|32|) can be expressed like 

zP2k = Z~'^{P2k+2 - 0-2k+2{k)P2k+l) + - Oi2ka2k+lP2k " (1 " |a2fe P)a2fc+1-P2*fc-1 > (38) 
^P2k+l = Ot2k+lP2k+2 - a2k+ia2k+2P2k+l + ^((1 " \a2k+lf)a2kP2k + (1 " |a2fc+lP)(l - |a2fc|^)-P2*fc-l)i 

(39) 

relations that can be obtained also with the classical Szego recurrence formulae and their reciprocals in T. 
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2.4 Projection operators and the Christoffel— Darboux kernel 

Here we discuss the CD kernel; i.e., the integral kernel of the quasi- orthogonal projection, according to the 
sesquilinear form (•, defined by the measure /i, to the space of OLPUC. 



Definition 8. We use the following notation 



^ ^ \Am], l = 2k + l. 

Notice that 

Associated with these spaces of truncated Laurent polynomials we consider the following related spaces, 
quasi- orthogonal complements, 

/<A:-Coo i<fc<oo 

Formally, we can express the following bi- quasi- orthogonality relations 

(AW,(A"^)^^)i- = 0, ((AW)^^AH)^ = 0, 

and the corresponding splittings 

A[^] = aW e (aW)^i = aW e {a^^Y\ 

induce the associated quasi- orthogonal projections 

^P:A[^]^AW, 4'):A[^]^AW. 

The reader should notice that we cannot properly talk of an orthogonal complement and an orthogonal 
projection if the measure is not positive and consequently we do not have an scalar product. If is a 
positive measure then (AI'])-*-! = (AW)-'-^ = (AW)-*- and both projections are truly orthogonal and coincide. 

Definition 9. The CD kernel is defined 

kV\z,z'):=Y,^['\z')^^^\-z). (41) 
fc=0 

As the CD kernel is expressed in terms of Laurent polynomials the definition makes sense as long as 
z,z' 7^ 0. This is the kernel of the integral representation of the projections 7r['\ VTg ^: 

Proposition 16. The integral representation 



holds. 



{4^f){z) = i K^Hz,z')f{z')d^^iz'), V/ G A[^], 



^In case that we have a positive measure fi then we can define the orthonormal Laurent polynomials (^''' = (hi) ^ipi'^ 
{h,)h^^^ so that kV\z,z') = Y.{':oK'^f\^')^f\-^) = Y!u=,h,^2\-')^f\-^) - EL=o^<''(^')^(--)- 
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Proof. It follows from the bi-orthogonality condition ([8]). □ 

This CD kernel has the reproducing property 
Proposition 17. The kernel K^^\z, z') fulfills 

(n, z')dfi{u). 

JT 

Proof. See Appendix Rl 

□ 

2.5 Associated Laurent polynomials 

In order to find CD formulae for the CD kernel just discussed we need the following definitions introdcing 
what we call associated Laurent polynomials 

Definition 10. In the I even case the associated Laurent polynomials are 



while for the I odd case they are 



The associated polynomials can be expressed in terms of the Laurent polynomials in difi^erent alternative 
manners 

Theorem 3. If is a positive measure the associated Laurent polynomials in the I even case are 
^fl,{z) = ^^l\z), ^f'_'^{z) = hi}Jl-'\z), 

<pll,iz) = z-^hi^^\z-^), = (42) 

= z-\aihi^f{z) + hi_,pf^t'\^)), = ^-Hphfi^) - c^i4~'\^)l (43) 

= /i/+i4'+'^(z) - hm+,^f{z), = ai.,^t'\^) + ^t'\^), (44) 

and in the I odd case 

^fl.iz) = ^f\z~'), 4;_'^{z) = hl\^f-'\z-^), (45) 

= z(a,^«(z) + p,Vf-')(z)), = z{hY^p^i^^P{z) - aih-_\^t'\^)) (46) 

= ^f+'\z) - ai+,^f\z), = hl\ai.,^t'\^) + K\ipt'\z). (47) 

Proof. 1. To prove (l42]l and (|45]) we proceed as follows. On the one hand, when I is even (j48l) implies 



.(0 



and on the other hand, due to the Hermitian property of the scalar product, it follows for that 
^^\z-^)z-^dp{z) = 0, J- = _1 + 1, . . . , 1. 
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Hence, zip^2+2 ^ ^[1/2,1/2] ^-^^d solves the same linear system of equations that tp2\z ^) S ^[1/2,1/2] 
does. Consequently both Laurent polynomials are proportional. The equality is obtained from the 
coefficients in the power z~2. In a similar way, from ([M]) we see that 

iz-')dfiiz) = 0, J = -1 + 1, . . . , 1 - 1, / z-2z-'^t-ii^'')Mz) = 1, 

this means that z^jg -i^(-z) S ^[1/2-1,1/2] has the same orthogonality relations and normalization con- 
dition that ^^(z^^) ^ ^[1/2-1,1/2]^ so they coincide. Analogously, in the odd case we obtain ()15]) . 

2. For (j43p and (|46p we argue in the following manner. Using orthogonality relations for (ff\z~^) and 
ip2 ^\z^^) we conclude that 

ip^l\z''^) e span{v9f\v9f""^^}, 
<^2~^^(^~^) ^ span{v?^'\v?^'"^^}, 

and identifying coefficients 

^'/^2,+2(^) = = ^t\z'^) = aiiffiz) + p;Vf"^^(2;) = aihiipf{z) + pfhi_iip^2'^\z), 

(1-1), X -{l-l)t -In 2 (0/ \ \ 

that concludes the proof of ()43p . (j46p follows similarly. 

3. Finally, we proceed now to prove (j44p and (|47p . For the even case we compute the following integral 

xW(^)4!V2(^)d^W = ix^Hz){x^'+^\z))^d^^{z)- 

( 9o,i+i \ 

X^Hz)x^Hz)^dKzm-' 

\9l-i,l+i/ 

{90,1+1 9i,l+i ■ ■ ■ 9l-i,l+i) — {90,1+1 9i,l-i ■ ■ ■ 91-1,1+1) 
(0 ... 0)^, 



which written componentwise reads 

£ z^^^l^mi^i^) = 0, J- = _!,..., 1 _ 1. (48) 

It also follows from the definition that +2(-^ ^ ^[i-i ((^2+2 ~ z^^^) g Ajj__-^ Fo^' 

the other associated Laurent polynomials, the orthogonality relations are 

z^^t-ii^Wiz) = 0, j = -1, . . . , 1 - 2, / z-2-'^t-imKz) = 1. (49) 

To get this result we proceed as before 

X^^z)^t-imKz) = ((fx^Hz)x^Hz)^df,iz)) (gW)-^e;_2 = e,-2, 

T \Jt J 

that is the matrix version of the orthogonality relations. 
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These orthogonality relations lead to 



Let us prove this statement. As iP2\2 ^ --i] t^i^n (/52 +2 ^ ^P^^i^'i'' ^^2^' ^ • • • ' '^2^^''}5 t)ut due 

to the orthogonality relations all the coefficients vanish except for the ones corresponding to ^2^^^ 
and (/?2^. Comparing the coefficients of z~^~^ and Z2 we get the system of equations 

1 = (5'2)/+\,/+ia/ + 0, = ai{S2)r+i,i^iai+i + {S2Xlk, 

from where we conclude ai = {82)1+1^1+1 a-nd hi = —ai+i{S2)i^i- Now, we notice that V'2-1^ ^ 
span{ip^\ ip^\ . . . , (^2 ^iid also that the orthogonality relations imply that 



(l-i) (/-I) , , (1^2) 

f2,-l = (^i-iV2 +dl-iip\ ■ 



(0 



.(0) ,.(1) 



^2,-1 -Lspanlv?^""^ , v?^^^ , . . . , <^2 }• 

Therefore, 'P2-1 ^ span{(/32 '^\^2 ^^'^ o^^y need to find the expression of the associated 
Laurent polynomials as a linear combination of these two Laurent polynomials. For that aim, we take 
the complex conjugate, multiply by ipf and (p^ and integrate to obtain 



(0) ,.(1) 



Q-1 



d, 



iff ^\z)(p^2-Wz)dfi{z) = f z 2(^^'_'^^(z)d^(z) + Q,_i ^ z2 ^(f':^_\>{z)dfi{z) = + ai_i, 



l~i = f ^1 '^\z)^pI_^^ {z)d^{z) = 6 Z2 'ip'^^ \\z)dn{z) = 1, 

so we conclude q_i = ai^i and = 1. For the odd case one proceeds in an analogous form. 
A different proof for the same formula can be found in Appendix [XI 

Finally we give determinantal expressions for these polynomials 
Proposition 18. The associated Laurent polynomials have the following determinantal expressions 



□ 



det qV] 



det 



50,0 50,1 • 


• 5o,«-i 




51,0 51,1 • 


• 5i,i-i 




5«-i,o 5z-i,i • 


• 9i-i,i-i 




5/+a,o 5z+a,i ■ 


■ 9l+a,l-l 





\ 



I > I. 



(50) 



and 



■1) 



det gV+^ 



det 



50,0 50,1 • • • 5o,«_a-i 
51,0 51,1 • • • 9i,i-a-i 



\9i,o 91 



( 90,0 
51,0 



9l,l-a-l 



det flW 



det 



50,1 
51,1 



5/-i,o 5/-i,i 



50, z_a+i ■ ■ ■ 9o,i 

51, «_a+i • • • 5l,i 

9l,l-a + l ■ ■ ■ 9l + a,l 



X^'\z)\ 

x'-'Hz) 



X^'\z) j 



9o,i-i 
91,1-1 

9l-i,l-i 



WHz))^ ix^'Hz))^ ... ix^'-'Hz))^ 



90,l + a \ 
9l,l + a 

9l-l,l + a 



{X^'^^\z))^ I 



I > I. 



(51) 



/ > 1. 



(52) 
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_(/) 



-1) 



l+l- 



det c/I'+i] 



■det 



/ 50,0 50,1 
51,0 51,1 

9l-a-l,0 9l-a-l,l 

9l-a + l,0 9l-a + l,l 



9L0 



91,1 



90,1 \ 
9i,i 

9}_a_-_}£ 

9l-a + l,l' 

9l,l 



I > 1. 



(53) 



Proof. See Appendix Rl 



□ 



2.6 The ChristofFel— Darboux formula 



To obtain CD formula in this context we need a number of preliminary lemmas. First we consider a version 
of the Aitken-Berg-Collar theorem |59j 

Lemma 1. The following ABC type formula 
is fulfilled. 

Proof. See the Appendix lAl 



The CD formula can be obtained using the previous expressions for the CD kernel. 
Lemma 2. For the CD kernel one has 

{z' - -z~')K^\z, z') = x^'\z)\g^'^)-^z'x^'\z') - -z-\^'\z)\g^'^)-\^'\z' ) 

= (X"l(^)n5f'l)-^5f''^'l - xf^'l(^)^)T[^'''k5['l)-VW(^') 
- xf'k^)n5f'l)-^T[''^'](5[^'''l(5W)-VW(^') - x^^-\z')). 

Proof. See Appendix Rl 

The reader can easily check 
Lemma 3. If I is an even number 

while for the I odd case he have 

Theorem 4. The following CD formula holds 



■1,1-2 = eoe;_2- 



K^'\z,z' 



(1 - z'z) 



zz / 1. 



□ 



□ 



(54) 
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Proof. The proof of (I54p rehes in Lemmas [2] and O Let us first study the I even case; with Lemma [2] and [3] 
we obtain a more exphcit expression for the CD kernel given by 

then using Definition [TU] for the associated Laurent polynomials and Proposition [3] we conclude our claim, 
that leads to ([Si]) . 

For the I odd case, reasoning again with Lemmas [2] and [3] we obtain the expression 

{-z-'-z')K^z,z') = -xW(-)^(5['l)-^5"'^'leo)e^_2(^7^'l)-^xW(-')- 
and recalling with Definition [10] we immediately get the claimed result. □ 

Recalling the different expressions for the associated Laurent polynomials in Theorem[3l one easily notice 
that 

Corollary 1. For a positive measure fi the CD kernel can be written in the following alternative forms. In 
the I even case it can be written as 



(l-z'z) 



iai^f\z) + pyi-'\z))^t'\^') - V?{z'){pf^^^\z) - aivt'\^)) 

(l-z'z) 



z{^f+'\z) - ai+,^f\z))^t'\^') - ^f\z')z{ai^,^t'\^) + 



(55) 
(56) 
(57) 



(1 - z'z) 

while in the I odd case it can be written as 

(1 z z) 

_ - Ml\-z){pf^^^\z') - ai^t'H^')) - io^i^?iz') + p^y^t'\^'))^t'\-^) 



z^?{z){ai^^^t'\^') + V>V\z')) - {^f^'\z') - ai+,^f\z'))zipt'\z) 



1 - z'z) 



(60) 



Formulae ()56p and ()59p where founqjin [25], however the other expressions, to the best of our knowledge, 
are new. 



3 Extended CMV ordering and orthogonal Laurent polynomials 

This section is devoted to an extension of the CMV ordering that allows to extend CMV results to more 
general situations. The main result is an extension of the CD formula allowing for projecting kernels to 
general spaces of Laurent polynomials, avoiding the CMV restriction on degrees. 

^However, the authors use the orthonormal sequence instead of the dual monic orthogonal sequences 
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3.1 Extending the CMV ordering 

Let us consider a vector n G Z^, n = (n_(_,n_), and the associated alternated sequence with n+ positive 
increasing powers of z fohowed by n_ negative decreasing powers of z,that is 

Xn{z) := ihz,..., z^+-\z-\z-\ z-^- , , ...y. 

The CMV case presented above corresponds to the particular choice n+ = ri- = 1. Given / > if is 

a non-negative power of z we say that a{l) = 1 and if Xfi^ is a negative power of z we define a{l) = 2. In 

addition, for any / > we will denote as the number of elements in the set {x^ ,0 < I' < 1} 

with a{l') = 1 {a{l') = 2). That is 

Ml) ■■= #{xP,a(0 = 1,0 <l'< I}, := #{x^p,a{l') = 2,0 < l' < I}. (61) 

We will use the notation 

z7 = {v^, := + |n| := n_|_ + n_. 

Additional expressions for ()6ip can be found using the Euclidean division, ([12). [16j). as we now explain. 
For any given / > there exists unique non-negative integers q{l) and r{l) such that 

I = q{l)\n\+r{l), < r(/) < n+, if a{l) = 1, 

I = q{l)\n\+n++r{l), < r(/) < n_, if a{l) = 2, 

so that 

^ (n = + + 1, a(0 = 1, ^ iq{l)n-, a{l) = 1, 

1(^(0 + 1)^+, a(0 = 2, \g(On_+r(0 + l, a{l) = 2, 

from where = / + 1. If {ek}'kLo is ^^^^ formal canonical basis of we consider {eaik)}'kLoy with 

o = 1, 2, defined as 

ei(z^+(/) - 1) := ei, e2{v-{l) - 1) := e;, 

these are new labels for {e^j^Q adapted to n. Given a non-negative integer / there exist a unique non- 
negative integer k and a number a G {1, 2} such that ea{k) = ei. This basis allows for a natural decomposition 
of Xn using positive and negative powers. In particular 

XnA^) ■=^ea{k)z^ , a = 1,2, Xn = Xn,l + Xh,2- 

k=0 

With those sequences we define the extended CMV moment matrix 
Definition 11. The extended moment matrix is the following semi-infinite matrix 



9n ■-- 



Xniz)Xniz)^df,iz). (62) 



T 



Notice that this moment matrix is a definite positive Hermitian matrix if is positive. The Gaussian 
factorization for the semi-infinite matrix (7^ is 

9n = {Sfi,l) ^Sfi^2, 
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where Sft^i is a normalized lower triangular matrix and S^^2 is an upper triangular matrix. The associated 
sequences of Laurent polynomials are 



where 



^n,l(2) 



V 



V 



As in the CMV case, the sets of Laurent polynomials {v^I\}'^q and {^^Ji\}'iZo bi-orthogonal with respect 
to the sesquilinear form (•, ■)_^, that is 



l,k = 0,l,.... 



(63) 



In addition if the measure fj, is positive we have that = hj^ ^^fi\ both families are proportional 
Laurent polynomials. In addition 



.(0 



/,A: = 0,1,... 



so in this case {v^'\};^o ^ orthogonal Laurent polynomials with ip^^^{z) G ^[u-{i),v+{l)-i\- The 

orthogonality relations (|64p read as follows 



(64) 



{V%z^)j^= f <^i';^(z)z-MM^) = 0, fc = -z._(/-l),...,z.+ (/-l)-l. 



.(0 



(65) 



In terms of the Szego polynomials we have 



Proposition 19. For a positive measure fi we have the following identifications between the extended CMV 
Laurent polynomials, the Szego polynomials and their reciprocals 



z-«4'\(z)=Pr(z), a{l) = 2. 



Proof. See Appendix [Al 



(66) 



□ 



The reader should notice that in the CMV case, n+ = n_ = 1, for I = 2k we have a{l) = 1, y-{l) = k 
and vj^[1) = A; + 1, and when I = 2k + 1 then a{l) = 2, V-{1) = = k + 1. 

3.2 Functions of the second kind 

Here we just give a brief description account of this extended case 

Definition 12. The partial second kind sequences with the extended ordering are given by 

Cri,l,l{z) ■■= {SZ^^yXH,l{z), Ch,1,2(.z) := (SZ^^)^ Xn,2{^) , CfiXl{z) ■= Sji,2XH,l{z), ChX2{z) ■■= SH,2Xn,2(.z). 

and the second kind sequences 

Cn,i{z) := (.S^DhHiz), Cn,2{z) := Sh,2X*h{z). 
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Generalized determinantal formulae can be obtained 
Proposition 20. The extended second kind functions have the following determinantal expressions for I > 1. 



/ 9n,0,0 9n,0,l 
9n,lfl 9n,l,l 



det 



9n,0,l \ 
9nA,l 



9n,l-l,0 9n,l-l,l ■ ■ ■ 9n,l-l,l 



f{0 f(0 

y n,2,0 n,2,l 



p(f) 



(67) 



and 



det 9^y 



■det 



^ 9nflfi 9n,0,l • • • 9n,0,l-l 
9n,l,0 9n,l,l • • ■ 9n,l,l-l 

\ 9n,l,0 9n,l,l • • • 9n,l,l-l 



n,2,l ) 



r(0 \ 

p(0 
^n,l,l 



r(0 , 



(68) 



where TV^^^. := Y.k>i9n,ikX^ and T'J^j := Y.k>i9n,jkXH ■ 

The same can be said about the relationship between the second kind functions, the Fourier series of the 
measures, and the integral representation, that can be found in 



Proposition 21. The partial second kind functions can he expressed as 



|fe|<oo |A;|<oo 



|fe|<oo 

and the second kind functions as 



and in 

Proposition 22. Assume a positive measure ji or a complex measure w{6)d6 with w a continuous function. 
Then, the second kind sequences can be written as the following Cauchy integrals 



|fe|<co 



(0 



nil) -1 1 "'^k2("):r^ 



\z\ > 1, 



\z\ < 1. 



3.3 Recursion relations 

As we already commented above the recursion relations among extended Laurent polynomials are more 
involved than in the CMV case. 
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Definition 13. Given n € we define the projections 



a = l,2, 



k=0 



and the shift matrices 



^n,a ■= ^ ea{k)ea{k + 1) 



A;=0 



T\n+ 



In the context of section [3] recursion relations can be obtained using the same technique. Our objective 
is to have an expression for the multiplication by z and by using the shift operators. 



Proposition 23. 1. The moment matrix g^i has the following symmetry 

J- n9n 9n ^ n • 



(70) 



2. The operator of multiplication by z has a |n| +3 diagonal band structure in the basis given by or 



Proof. See Apendix[ 



□ 



Now we introduce the following associate integers. For a = 1, 2 we shall call Z+a and l-a to the smallest 
(largest) integer /' that verifies I' > I {I' < I) and ai(/') = a. For instance, in the previous case with 

= = 1, if I is an even number = /_i = /, /+2 = 1 + 1, and /_2 = / — 1, in the case that / is an odd 
number then /+2 = 1-2 = I, l+i = 1 + 1, and l-i = 1 — 1. This structure leads to the following recursion laws 
for k >1 and < Z < |n| — 1 (the A; = case corresponds to the truncated recurrence relations). 



z(p, 



(|n|A:+0 



n,l 



((|n|fc+«+l)+i) 



((ln|fc+i+l)+i-l) 



+ • • • + 4 



{{\n\k+l-l)^2) 



where mfi{l) is the number of non vanishing terms in the recurrence formulae. The coefficients j{k) 
are again labeled with the index j that accounts for the mf({l) non vanishing coefficients for each I; as 
there are only \n\ "different" recursion relations (the equivalent to the recurrences for the odd and even 
polynomials) then / = 0, 1, . . . |n| — 1. The connection with the elements of the Jacobi operator is the following 



Jl 



J, 



n,\n\k+l,{\n\k+l). 



The reader can check that m^(/) = {\n\k + I + l)+i — {\n\k + I — l)-2 + 1- 



Due to the fact that {\n\k + I + l)+i < \fi\k + / + n_ + 1 and that {\n\k + I — 1)_2 ^ l^^l^ + / — 1 — n_|_ then 
— 1^1 + 3 that agrees with the structure of \n\ + 3 diagonals. Furthermore it is possible to calculate 
mfi{l) more explicitly and show that actually it does not depend on k. We can see that 



+ / + l)+i 



\n\k + l + l, <l <n+ -2, 
\n\{k + 1), n+ - 1 < / < |n| 



and consequently 



\n\k + l- 1) 



\n\k — 1, < I < n^, 

|n|A; + /-l n+ + 1 < / < |n| - 1, 



mfi{l) 



1 + 3 l = 0,...,n+-2, 

|n| + 2 / = n+ — 1, n+, 
\n\ — I + 2 / = n+ + 1, . . . , |nj 
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so in fact ra^{l) < \n\ + 2. 

The expressions for the coefficients J^p{k) with / = 0, . . . , n+ are 



iSfi})i\n\k+l+i)+u{\n\k+l+i)+i-p P = 0, . . . , {\n\k + / + l)+i - {\n\k + /) - 1 



and for I = n+, . . . , [n| — 1, the expressions are 

jl ^ f(5'^,l)ln|fc+«,(ln|fc+0-i('S'^,i)(|n|fc+/)+i,(|n|fc+0+i-P P = 0, • • • , (|n|A; + - {\n\k + /) 
"''^ \ {hfi)\H\k+l{S^^2){\n\k+l-i)-2,{\ft\k+l)+i-p p={\n\k + l + l)+i - (|n| /c + / + 1) . . . , - 1. 

The particular case of n_ = n+ = 1 gives the 5 diagonal CMV matrix with only four non-vanishing 
coefficients. As a consequence, the standard CMV case has the shortest possible recurrence formula. 

3.4 The Christoffel— Darboux kernel 

We discuss the CD kernel for this extended case 

Definition 14. For each non-negative integer I we define the set of truncated Laurent polynomial subspace 
as the following span 

:=C{xf ,...,xr'^} = A[..a-i),.,a-i)-i]. 
As before, we define quasi- orthogonal subspaces 

(aS)^^:={ E c,^^&c,ec}, (Aff)^^:={ E ^^^'U. € c}, 

«<fc<oo /<fc<oo 

so that the following bi- quasi- orthogonality relations are satisfied 

(Af , (A?)^Oi? = 0, ((Af Af )^ = 0, 

and the corresponding splittings 

A[^]=A|!l©(AW)-^=AU©(A|!l)-^ 
that induce the associated projections 

vr^'^, : Ar^i ^ A^, Tr^!\ : Ar^, ^ a'^, 

n,l [ooj n ' n,2 looj n ' 

hold. In the positive definite case this extended version allows for the projection in more general spaces of 
truncated Laurent polynomials. Recall that for the CMV case the space of truncated polynomials given in 
(j40p includes only a very particular class of these truncations, excluding the majority of cases. The intro- 
duction of the extended ordering allows to include in the discussion all the possible situations of truncation. 
In fact the space A[p g] can be achieved in a number of ways, and always by the choice n_|_ = g + 1, n_ = p. 

Definition 15. The CD kernel is 

i^W(.,z'):=E4';l(-')4';^(--), m) 

A:=0 

and, whenever the measure fj, is positive definite, we have the equivalent expressions 

/-I i-i i-i 



fe=0 k=0 k=0 
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Proceeding as in the CMV ordering we conclude the fohowing results 
Proposition 24. 1. This is the integral kernel of the integral representation of the projections tt^\, Tr^\ 

(4?i/)(^') = / K^^{z,z')f{z)di,{z), V/ G A[^], 

Jt 

2. This CD kernel K^3{z,z') has the reproducing property 

kI\z,z') = I k'§{z,u)k'§{u,z')Mu). 

JT 

3. The following version of the ABC theorem holds 

Kf{z,z')=^^{z)\9^l)~^^^{z'). 

4- The CD kernel has the following expression 

{z' - -z-')KiHz,z') = xff (.)t(5f )-^.'x!!^(^') - z-'^{z)Hg^)-'xS{z') 

The integers l±a can be used to calculate the T blocks in this case. The reader can check 
Proposition 25. The formula for T^j is the following 

Thus, the expression of the CD kernel is 

- z')K^^{z,z') = [:^\zy - xi(.)t(4')-^4'- ')e/..-^4-i)-.(^f 

that suggests the definition of the following associated polynomials 
Definition 16. The associated Laurent polynomials are defined by 

r^n,l,+a X^ \yn,l+a,0 yn,l+a,l "' yn,l+a,l-^) Kiln ' ■> y^n,l-a ■~ ^l-ay^n > ' 

JO ._Al+a)_(-^ - - ._pT // [i+l]wUt^[;+l] 

where a = 1, 2. 

It is easy to see that 4?i,+a(/) = ' 4' 2,+a(0 = ^^■^)^^^%^^%,-a{l) = ^-^^^ and 9'^^^2,-a(0 = ^'V 
Theorem 5. For the associated Laurent polynomials f^^-a have two alternative expressions 



(72) 
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1. The reciprocal type form (valid for positive definite cases) 
when a{l) = 1 and 



(73) 



(74) 



for a{l) = 2. 

2. The linear combination form (valid for quasi- definite cases) 

4'!i,-a = (^2-')/-.,^-.4V^ + (^2-')^_.,/-.+i4V^'^ + • • • + is:^')i-a^% (77) 

= {Sr)i_i-A:2^ + (^O^^.+i./^.^^""'^ + • • • + (^i)'.'-.4'2- (78) 

Proof. 1. Let us suppose that a{l) = 1. In that case we have ip^^^ ^-^ = 93^'^^ and ip^^-^ ^2 ^ ^[[/-(O+i.i'+W-s]- 

Consequently, (')+V!J?i,+2 ^ A[.+ (/)-i,._{0] and a(/) = 1, = u+{l + U-{1) = U-{1- 

1). For the dual polynomials -1 = '^'h,2 -2 ^ ^li^-CO.i^+CO-i]' ^^^nce z''- (')+V^'^2 -2 ^ 

-A[j/+(z)-i, ;/_(/)]• Using (j79|) we conclude that the following orthogonality relations hold true 

^._a)-.+ (0+2^W ^^^(^)^-fcd;,(^) = 0, A; = -^.+ (/ - 1) + 1, ...,!._(/ - 1), 

^._(0-.+ (0+i<^W^_^(^)^-fcd/,(z) =0, k = -u+{l-l) + l,...,u_{l-l), 



IT 

/T ' ' 

and we get the result. 

Now let us suppose that a{l) = 2. In this case we have 9?^'^ ^2 ~ v'ni ^"^^ '^^Hi+i ^ 
Consequently, z^-^^^~'^+^^^Lp^^^ € ^[v+{i)-i,v_{i)]- Now, as a{l) = 2, we have = — 1) and 

i/_(/) = u^{l - 1) + 1. For the dual polynomials ^^^^2-2 = '^i',2 ^""^ '^i',2-1 ^ A[i^-{0,i'+{0~i] ' ^° 
^^_(0-^+(/)+i(^W^ _^ G A[^^(;)_i^^_(,)]. Now using again (HI we get 



^^^(^)^-fcd^(^) = 0, A; = -z.+ (/ - 1) + 1, . . . , z._(/ - 1), 

= 0, A; = -z.+ (/ - 1) + 1, . . . , z._(; - 1), 



2. For ^^^i direct computation gives 



4'!i,+a(^)xi(^)^d^(z) = iixl^-Hz) - {g^ua,0 9nUa,i ■ ■ ■ 9HU.,i-i) (^f )-'xi (^))xS (^)^dM(^) 



{9n,l+a,0 9n,l+aA ' ' ' 9n,l+a,l-l) ~ {9n,l+a,0 9n,l+a,l ' ' ' 9n,l+a,l-l) 

0, 
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and for wl'^ we have 

IT JT 

so that we get orthogonahty relations for the associated polynomials 

+ai^)z'''Mz) =0, k=-l^.{l-l),..., U+{1 - 1) - 1, 

(fc) 



(79) 



Therefore, 



Xh {z)(Pn,2~a{z)<ifJ'{z) =^fc,/-a' k = 0,l,...,l. 



(0 

I.e., Lp-^\ , 



X]jt=l ^j'V^i fo'^ ^ set of coefficients {A^-^}. Comparing the powers of z that appear in 
the subsequence {x^-'^}i<j<i+a ™ both sides of the equation, the following linear system of equations 
is obtained 



/I 

(5'l)«+„,i+„-2 (5'l)i+„-l,i+„-2 1 
\ (5'i)/+a,/ {Sl)l+,-l,l 










o\ 







('S'l)i+2,i {Si)l+i^l 1/ 



t + a — i 
'+a— 3 



/1\ 







calling Ai the coefficient matrix, the solution can be written as 



/ 4^ \ 



\ A'" / 



Prom the structure of A4 we conclude that 



/O 


1 
VI 



1\ 

1 



0/ 



M 



/O 




1 

VI 



1\ 

1 


0/ 



V(M-i)/+„-V 
/ 1 



(Si) 



1+1,1 





1 







V iSl)l + a,l ('S'l)i + a,/ + l 



1 



From the triangular structure of Si we deduce that {A4' = {S^ for i, j = 0, 1, . . . , l+a — I 



and consequently {M. = {^A 



)l+a-l,l + a—l-j 



l+a,l+a—j 



for i,j = 0, 1,...,/ 



L which 



proves ()75p . The expression for (|78p is obtained using a similar technique. Using again ()79p we conclude 

that 4'j2,-a e span{4'7\4'7+^\...,4')2}; i.e., 4\„„ = E5=/_. V'S/^l Bi-orthogonality and 
normalization properties imply 



B 



(0 



j = l- 



-a; • • • ) ""J 



that proves ()78p . The other two equations are obtained using the same idea. 



□ 
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The polynomials that appear in the CD formula are now clearly identified as 
and consequently the final result is the following 

Theorem 6. The CD formula for the extended ordering is the following 

We get the followings corollaries when we have a positive Borel measure jj, 
Corollary 2. Given a positive measure jj,, the CD kernel can be expressed using 

in the case a{l) = a{l — 1) = 1, 





(1 - z'z) 




1) = 1, 














(1 - z'z) 




- 1) = 2, 








-i)_^'-+(0---(O^W^, 


{z'-^)-z4;^\-z) 




(1 - z'z) 


•> 










'-ii)^f^{z'-^)z'^+^i)-^- 


(o+vfei)(,-i) 




(1 - z'z) 


5 



KlKz,z') 

in the case ail) = a{l — 1) = 2. 
and 

Corollary 3. The CD formula for a positive Borel measure /x can be expressed in terms on the Szego 
polynomials as 

K^^iz, z') = /,-_i^^a(0+^+(0-2^/a(0-i-^-(0 ^K^-^)^r-i(^Oj-^i^(^O^r-i(^-^) _ ^33^ 

4 Associated 2D Toda type hierarchies 

Here we analyze the link between the previous constructions on OLPUC and integrable systems of Toda type. 
Our driving idea is the presence of the Borel-Gauss factorization problem in the theoretical construction of 
both, OLPUC and Toda. 
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4.1 2D Toda flows 



We will consider a set of complex deformation parametera^i t = {tij,t2j}jeN and with it two semi-infinite 
matrices that we will define now0 

Definition 17. 1. The deformation matrices are defined as follows 



Wiflit) := exp ( J2 ' ^2,o(t) := exp ( ^ t2j{T^yj , (84) 

i=i i=i 

2. For each t we will consider the matrix g[t) 

g{t) ■.= Wi,^{t)g{W2fi{t))-\ 

3. and the corresponding time dependant Gauss-Borel factorizatioi^^ 

g{t) := Wifl{t)g{W2,oit)r\ g{t) = {S,it))-'S2it). 

As we show now the deformed moment matrix is a moment matrix of a deformed measure. 

Proposition 26. The "deformed" moment matrix can be understood as a moment matrix for a "deformed" 
(that is, a "time" dependant) measure given by 

oo 

d/x(t, z) := exp (^''^^tijZ^ — t2jZ~^^dfj,{z). (85) 
Proof. First expanding the exponentials in (184p we obtain 



k=0 

then using the definition of g and g{t) we get the desired result 



oo „ oo oo 

k,l=0 •^'^ k=0 1=0 

„ oo 

= f Xiz)xiz)^ (i^P (^^{hjz^ -t2jz-^)^dfi{z). 
From this result we conclude at least for absolutely continuous measures 

oo 

F^{t) = exp (^^tyz^ - t2jZ~^^F^{z), 



□ 



^''in the framework of the theory on integrable systems these parameters are understood as an infinite set of times, being the 
independent variables in an associated nonhnear hierarchy of partial differential-difference equations 

^^We shall drop the subindex n from g and T as the definitions are valid for any value of n. It will be supposed that a 
particular n is chosen and the whole iJ3]will be built using that n. 

^■^For the sake of notation simplicity in some situations we drop the time dependence of Si, S2 and they will not denote the 
factors within the Gauss-Borel factorization of the initial condition but for the "deformed" one. Consequently, in this section 
Si , S2 will always depend on "time" parameters. 
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from where we deduce that the radii defining the annulus of convergence is time independent; i.e., R±{t) = 

R± 

Given a positive definite initial measure fj, in order to ensure that the evolved measure ^{t) is also definite 
positive for all times is enough to request to the exponential to be real; i.e., setting t2j = —tij, so that 

oo oo 

exp '^{hjZ^ + tijz~^)^ = exp ^ 2 Re{tijZ^)^ . 

j=0 j=0 

4.2 Integrable Toda equations 

Definition 18. Associated with the deformed Gauss-Borel factorization we consider 

1. Wave semi-infinite matrices 

Wi{t) := Si{t)Wi,o{t), W2{t) := S2{t)W2,o{t). (86) 

2. Partial wave and partial adjoint (denoted the adjoint by *) wave semi-infinite vector function^^, 



Wi{t)xi{z), n,i{z,t) 

Wiit)x*2iz), n,2{z,t) 

{Wi{t)-^)^xl{z), ^2,l(^,t) 

{Wi{t)-')^X2{Z), ^2,2{^,t) 



(87) 



^i,i{z,t) : 
^i,2(-z,t) : 
n,iiz,t): 

and wave and adjoint wave functions 

^i{z,t) := Wiit)x{z) = (*i,i + ^i,2){z,t), %iz,t) := {W2{t)-')^x{z) = (^2,1 + ^2,2)(^' 

^Uz,t) := {W,it)-')^x*{z) = i^l, + n,2){z,t), Mz,t) := W2{t)x*{z) = {^2,1 + ^2,2){z,t). 



{W2{t)-')^Xi{z), 

{W2{t)-')^X*2{Z), 

W2{t)xl{z), 

W2{t)X2{z), 



3. Lax semi-infinite matrices 

Li{t) := Si{t)rSi{t)-\ 

4- Zakharov-Shabat semi-infinite matrices 



L2{t) ■=S2{t)T' S2{tr'. 



B2,, := {Liy 



(90) 



where the subindex + indicates the projection in the upper triangular matrices while the subindex 
the projection in the strictly lower triangular matrices. 

Theorem 7. For j,j' = 1, 2, . . . ,the following differential relations hold 

1. Auxiliary linear systems for the wave matrices 



dWi 
dtij 



dWi 
dt2j 



B2M1, 



dt 



B1JW2, 



dt 



2j 



2. Linear systems for the wave and adjoint wave semi-infinite functions 



dtij 
dtij 



dt2j 
dt2j 



on 

dtij 



B2M2. 



(91) 



2'-? ^' dtij 



2j 

d^2 



dt 



2j 



B2,j^2- 



(92) 



Also called Baker, or Baker-Akhiezer, functions 
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Li], 




dBi,f 


dtij' 






dB2,j' 


dt2j' 


dt2j 


dBi,j 


dB2,f 


dt2j' 


dtij 



3. Lax equations 
4- Zakharov-Shabat equations 

+ [Bij,Si,,.] = 0, (94) 
+ [52,^,^2 J'] = 0, (95) 
+ [Bij,S2jv] =0. (96) 

Proof. The proof can be made using the same idea used in [T2], so we do not repeat them here again. □ 

From the definition it is clear that the wave functions are associated to the OLPUC for the evolved 
measure 

Proposition 27. The wave functions are linked to the OLPUC and the Fourier series of the measure trough 

^t\z,t) = ^t\z^t)e^7^^'^r^\ 
(M/*)W(.,t) = #(z,t)e-Sr^i*^^-\ 
mt\z,t) = 27r4")(z-i,t)z-iF^(,)(z)e-^r^i*^^-^^' 
= 2v^4")(z-^^)z-lF^(z)e-^r^^*^^-~^ 
*(")(.,t) = 2vr4'^)(z-\t)z-iF,(,)(z-i)e^r-*^^-^ 
= 27r4")(z-i,t)z-iF^(z-i)e^r=i*^^^"' . 
Moreover, the wave functions are eig en- functions of the Lax and adjoint Lax matrices 

Li^i = z^i, Ll^2 = z'^h 

LI^I = Z^I, L2^2 = Z^2- 

4.3 CMV matrices and the Toeplitz lattice 

For the CMV ordering of the Laurent basis, Lax equations ()93p can be we written as a nonlinear dynamical 
system that is a version in the CMV context of the Toeplitz lattice discussed by Mark Adler and Pierre van 
Moerbeke [8] 

Proposition 28. For the case n = (1, 1) Lax equations ()93p have as a consequence the following nonlinear 
dynamical system for the Verblunsky coefficients 

- .,(1) n _ .,(1)^,(2)^ da^ _ __^,(2) _ Ji)^,{2)^^ 

(99) 



with k = 1,2, . . . . 
Proof. See Appendix Rl 



□ 
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If the initial measure /i is positive definite then the sequences {a^^''} and {a^^''} are identical in t = 0. 
Furthermore, if we set t2i = —tu, the evolved measure is always real and there is only one family of 
time-dependent functions. That is the reduction studied by L. B. Golinskii [36] in the context of Schur 
flows. 

We have obtained a CMV version of the Toeplitz lattice, but for any n the integrable hierarchy obtained 
is always equivalent to the one discussed in [8]. The reason for this fact relies in the observation that for 
any positive Borel measure /i and for any n, there is a bijection between the set of OLPUC {f/^^'^} and the 
set of OPUC {Pi}. All the coefficients of any Pi are determined in terms of the set of reflection coefficients 
{ai}, so the time evolution for the Szego polynomials under Toda-type flows is determined by the evolution 
of the reflection coefficients. As the measure evolution does not depend on n, it is natural to always obtain 
the very same evolution for the family {ai} under the Toda flows. We conjecture that a similar result holds 
for the quasi-definite case. 

4.4 Discrete flows 

We now consider discrete flows associated to the moment matrix. Given two integers si, S2 and s := (si, S2) 
is then possible to make a new deformation of the moment matrix that depends on s. 

Definition 19. We introduce for each s and the deformed moment matrix g{s) and its deformed Gauss- 
Borel factorization 

g{s) := Di,o{s)g{D2,o{s))-\ g{s) = S{\s)S2{s), 

where Di^q{s), D2fl{s) are discrete deformation operators to be determined later on. 

We consider the operator Ti responsible of the shift si si + 1 and T2 corresponding to the shift 
S2 ^ S2 + 1. Let us suppose that matrices qi,q2 exist and satisfy 

Tipi.o) = qiDi,o, Ti{D2,o) = ^2,0, 

T2{Difi) = Difi, 72(1)2,0) = 92-02,0, 

then we define 

61 := Si{s)qiSi{sy'^ , 82 := S2{s)q2S2{sy^ . 

If 81,82 can be LU factorized, then there exist semi-infinite matrices (52,+, (52,- such that 

81 = (5j;^i(5i,+, 82^ = (52^i(52,+- 

In this case we introduce 

wi := (5i,+ , UJ2 := 82-. 

Proposition 29. The operators Ti,T2 and the matrices Si{s), S2{s) satisfy the following equations 

Ti{Si{s)){Si{s)r' = Ti{S2{s)){S2{s))-' = 

T2{Si{s)){Si{s))'' = 82,- T2{S2{s)){S2{s))-' = 82,+ 

Proof. First using Ti and T2 on the factorization Di QgD2Q = S^^S2 we obtain 

TiiD,,o)gTi{D2^l) = Ti{Sy)Ti{S2) {Ti{Si)Sy)-^Ti{S2)S2' = 8^ 

T2{D,,o)9T2{D2^l) = T2{Sy)T2{S2) ^ {T2{Si)Sy)-'T2{S2)S2' = 82', 

then using the factorization for 81 and 82^ and its uniqueness we can identify the upper and lower triangular 
parts and prove the claimed result. □ 
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It is also possible to define wave matrices Wi and W2 in this discrete context, 

Wi := SiDi^o W2 := S2D2,o. 

To ensure the consistency between both continuous and discrete flows we only need to replace -Di.o 1— > 
DifiWi^o and 1)2,0 '-^ -D2,oW^2,o- The next results are valid in case continuous evolution is also present or 
not, for the proof one only needs a slight modification of the one in fl2j. 



(100) 



Theorem 8. 1. The next linear system for Wi and W2 is satisfied 

Ta{Wa') = iOaWa' a,a' = 1,2. 

2. The discrete versions of the Lax equations are the following 

Ta{La') = UJaLa'UJ'^ a, a' = 1,2. 

3. The compatibility equations for the discrete flows in the linear system (|100p are 

Ti{uJ2)uJi = T2{uJi)uJ2 

if there are also continuous deformation parameters the mixed compatibility equations are 

duja 



(101) 



(102) 



Ta{B2j) 



duJa 



dt 



2j 



uj^'^ + UaBiju^^ a = 1,2 j = l,2,... 
uj~^ + ujaB2jUj~^ a = 1,2 j = l,2,... 



(103) 



Now we give examples of some discrete flows operators. Let be {Ai(j)}jgz and {\2{j)}jei two complex 
sequences with \i{j),\2{j) € B, then 



D 



(1) ._ 
1,0 •- 



nH(T-Ai(-i)i 



^-1 



ni > 
ni = 
Til < 



2,0) 



-urn 



[n"4(TT-A2(-j)ii)- 



n2 > 
n2 = 
n2 < 



the evolution of the measure is then 

d^^{z,s) = V^^\z,si){v':^'^)-\z,S2)Mz) 

where 

'n^io(^-AiO')) ^i>o 

< 1 = {V^^^)~^{z,s) 

ni:ii(z-Ai(-j))-i .i<o 



V^^\z,s) 



[iqUz-^ - \2{j)) 



1 



[n%{z-^-\2{-j)) 



-1 



S2 > 
S2 = 
S2 < 



in that case 



„(i) 



T' - A2(S2 + 1)II 
L2 - \2{S2 + 1)1 



T- Ai(si + l)I 
5\'' = Li-\i{si + l)l 

The evolution of the wave functions is associated to the evolved Laurent polynomials 

^i(z, s) = W^{s)x{z) = Si{s)D,^o{s)x{z) = <^i{z, s)vf\z, s), 
n{^,s) = {W2{s)-^)^x{z) = {S2{s)-^)\D2fl{s)-^)h{z) = <^2iz,s){V^^Y\z,s), 
where ^i{z,t) and ^2(2, t) are the Laurent polynomials associated to the evolved measure. 
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Lemma 4. We have the following structure for the matrices u}i,u}2 

wi = a;i,o + H h t^i,n_+iA"-+-^ 

UJ2 = W2,o + ^2,1^^ + ■■■+ a;i,„++i(AT)"++^ 



Pifl + Pi,iA^ H h /9i,„_+i(A 



P2,0 + P2,lA H h pi,n++lA' 



TN,n_+l 



n++l 



for some semi-infinite matrices 



ujij = diag(a;ij(0),tJij(l), . . . 

uj2,j = diag{uj2,j{0),uj2j{l), ...[ 

Pi,j = diag(pij(0),pij(l), . . . ) 

P2,j = diag(p2,j(0),/)2,i(l),---) 

Proof. Immediate from (jl04p . 

Defining 

71(2;, s) := z- Xi{si + 1) 



j = 0,...,n_ + l 
j = 0, . . . , n+ + 1 
j = 0,...,n_ + l 
j = 0, . . . , n+ + 1 



□ 



72(2, s) := 2; - A2(s2 + 1) 



the previous Lemma allows us to compute the action of the operators Ti and T2 on the OLPUC ipf^ {z, s), Lp^2 {^^ 
Proposition 30. The following equations hold 

(Ti(^«)7i = ^i,o(/)^S'^ + + • • • + a;i,„_+i(/)v^l'+"-+') 

= pimTwf) + PiAi){TwV^) + ■■■ + pi,„_+i(0(riv.('-"-^)) 

= [p2mT2vf) + P2,imT2vt'^) + ■■■+ Pl,n^+l{l){T24^''^^'^))l2 

Proof. The first two equations come from (llOOp and Lemma |4l For the last two equations we use that 

^lU{W~,y) = {W-,y a,a' = 1,2. 

□ 

Another possible option that preserves the reality of the measure is using pairs of conjugate transforms 
as follows 



D 



(2) 
1,0 



2,0/ 





Ai(j)I)(TT 


-Ai(i)I) 




> 








Si 


= 


ni^i(T- 


Ai(-j)I)-^ 


(TT-Ai(-i)l 




< 




-A2(i)I)(T 


-A2(j)I) 


S2 


> 








S2 


= 


ni^i(TT 


-A2(-j)I)" 


i(T-A2(-i)] 


y S2 


< 



the evolution of the measure is then 



d/i(z, s) = Vf\z, s^){V'^>r\z, S2)Mz) 
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where 



n,^ol^-Ai(i)l' si>o 

1 si = 

tni^i|z-Ai(-j)|-2 si<o 



tt|s2| 1-1 



S2 > 
S2 = 
S2 < 



(104) 



in that case 

(2) 
a(2) 



(T-Ai(si + l)I)(T^-Ai(si + l)I) 



(2) 
,(2) 



(TT-A2(s2 + 1)I)(T-A2(s2 + 1) 



(Li - Ai(si + l)l){L^' - Ai(si + 1)1) 5),> = (L2 - A2(s2 + l)I)(i^' - A2(s2 + 1)1) 



(105) 



Observe that these discrete flows leads to extended Geronimus transformations [46j . When the sequences 
{-^i(j)}) {'^2(j)} are constant and thus (?i,(?2 are invariant under the action of Ti and T2 we can make an 
interpretation in terms of Darboux transformations in the context of |7]. In that case what we obtain is 



5i = 



1,+ 



T252^ = T2iW2)q2^T2iW^^) = 0J25^^0J^^ = 52,-52"i52,+52"i = 52,+52"i 



that is a change in the LU factorization into UL. 
4.5 r-functions 

As is well known r- functions is an essential ingredient of the theory of integrable systems. Not only for 
the use of Hirota of these functions in the construction of soliton solutions [38j but also for its relevant 
geometrical insight [26]-l28], also the bilinear equations discussed in the mentioned papers are fundamental 
in the construction of solutions. The determinantal expressions for the OLPUC, and the associated Laurent 
polynomials and the corresponding second kind functions lead to a r-function representation of these objects. 
For that aim one consider the action of coherent shifts in the time variables, the so called Miwa shifts 

Definition 20. 1. The Miwa shifts are the following time shifts 



2. And the Miwa dual shifts 



t t ± [w]i := tij l-T- tij ± t2j ^ t2j, 



t t± [w]2 := tij I-)- tij, t2j ^ t2j ± —. 



A very important property of this Miwa shifts is the form in what they act on the deformed measure 
Proposition 31. The evolved measure ^{z,t) has the following behaviour 



yL{z,t±[w-\) = {l--) fi{z,t), 
\ wJ 
/ ^\ ±1 

^i{z,t±[w\2) = [1 - —j fJ'{z,t) 



\z\ < \w\, 
\z\ > \w\. 

Proof. Using the series expansion of the logarithm, the evolution factors change under Miwa time shifts like 



exp 



00 00 ^ 

'^{ti,z^ - t2,z-^)) ^ eM^Hhj ± —-)z' - t2jZ~^)) 



1 ) expC}2{tijZ^ -t2jz~^)), \z\ < 



exp (^ (tij z-' -t2jZ ■^)) iH> exp (y^ (tij z-^ 



{t2,±-)z-^))=[l-- 



W 



±1 



exp(^(tlj2:-'' - t2j^; ■')), \z\> 



□ 
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We introduce the main and associated r-functions as determinants 
Definition 21. The r-function is 

r(0)(t):=l, T«(t) :=det5W(t), 

while the associated r-functions are 



rtaii) ■■= (-1) 



l+l-a 



det 



/ 50,0 50,1 
51,0 51,1 



5o,«_a-i ] 90,l^a+l 

gi,l-a-l ! 5l,«_a + l 



I 



\5i-i,o 5i-i,i ■■■ 5z~i,Z-a-i]5/-i,i_„+i 



/ 50,0 
51,0 



4iait) ■■= (-ir'-" det 



50,1 
51,1 



9o,i \ 
9i,i 



gi-a + 1,0 9l-a + l,l ■ ■ ■ 9l-a + l,l 

V 5/,o 5i,i ••• 9i,i J 



1 = 1,2,..., 



90,1 \ 
91,1 

9l-i,l J 



I = \n\, \n\ + 1, . 



/ = |n|, |n| + 1, . . . , 



f 5o,o 5o,i • ■ ■ 5o,«-2 
51,0 51,1 ■ ■ ■ 5i,«-2 



^5«-i,o 5z-i,i ■ • ■ 91-1,1-2 



I 50,0 
51,0 



Ti«,(t) :=det 



50,1 
51,1 



5/-2,0 5«-2,l 



50, (Z-l) + a \ 

51, (z-i) 

5;-i,a-i)+„ j 

90,1-1 \ 
9i,i-i 

9l-2,l-l 



I = \n\, \n\ + 1, , 



/ = \n\, \n\ + 1, . . 



V5(Z-i)+„,o 5a-i)+„i ••• 9(i-i)+a,i-i ) 
To find expressions for the OLPUC in terms of these r-functions we need the following 

of the matrix g{t), then for j € Z+ \ {0, n+j 



row 



Lemma 5. Let he rj{t) the j-th 

rj{t - [z-'li 

rj{t+[z]2) 



(t) - z a{j) = 1, 

a{j) = 2, 

^(t) - zr(j_i)_j(t) a{j) = l, 
'jit) - zr(j+i)+^{t) a{j) = 2. 



is satisfied. For j = or j = one has 

roit - [z-\) = roit) - z-^ri^,{t) 
roit + [z]2) = ro{t) - zro^^{t) 
rn+{t - [z-\) = rn^t) - z-Vo(t) 

rn+it + [zh) = Tn+it) - 

Proof. Immediate from Proposition [3TJ □ 
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Let us recall the skew multi-linear character of determinants and the consequent formulation in terms 
of wedge products of covectors. Observe that 

Lemma 6. Given a set of covectors {ri, . . . ,r„} it can he shown that 

n n+l 

l\ {zT, - r,+i) = n A r2 A • • • A f,- A • • • A r„+i, (106) 

where the notation fj means that we have erased the covector rj in the wedge product ri A • • • A r„,-|_i. 
Proof. It can be done directly by induction. □ 

These two lemmas are the key property to characterize deformed OLPUC using r-functions 
Theorem 9. Given I > \n\, one has the following r-function representation of the OLPUG 



for the other "+ "associated polynomials we have 



(107) 



(0 

and finally, the remaining "— " polynomials can be written as follows 

.(0 



(108) 



(109) 



Proof. Let us prove p07|) . If a{l) = 1 we can use Lemma [U] with ri = r^i_i'j_^ and r„, = r(^_]^)_^ to expand 

^.+(o-v(0(i _ [^-1]^) = ^-+(o-1(m(/+i) + (-i)'+('-i)-z-im(;+;j_^, + • • • + (-i)'+('-i)-^z-'Mf;+;)_^,) 

= ^(-l)'+^Mf^)x(^-)(.) 

j=0 

= det{g^Ht))^?iz,t) 
= r^^\t)^?iz,t). 

If a{l) = 2 the same procedure works with ri = r(;_i)_^ and r„ = r(;_i)_2- Now the expansion is 

a)^(0(i + = (0(a4'+i) + (-i)'+('-i)-2^mj;+;j_^^ + • • • + (-i)'+('-i)-z'Mf;+;j_^,) 

= ^(-l)'+^Mf^O-)(z) 

j=0 

= det(5W(i))rf^(^,i) 
= r«(t)v.«(z,t). 
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The proof of expressions in (llOSp can be performed using the same technique, expanding the right hand 
side. If a{l) = 1 we have to consider Lema[6]with the rows ri = r(^i_i'j_^ and r„ = r(^i_i^_,^. In the case 
a{l) = 2 the rows n = r^j^j^-j ^ and r„ = r(;_;^)_^ should be considered. 

To conclude we prove (|1U9|) . repeating previous arguments with adequate r- functions, that is, in the case 
a{l) = 1 



,+ 



j=0 



and in the case a{l) = 2 then 



□ 



To obtain the r-function representation of the dual Laurent polynomials (p2 ^ and their associated ones 
we would proceed using again Lemma [5] and Lemma [6] interchanging the role of rows and columns, that 
would lead to the final expression 

Theorem 10. For any I > |n| the dual Laurent polynomials 992 have the following expressions in terms of 
T -functions 

the "+" labeled associated polynomials can be written as 



to conclude, the "— " labeled polynomials have the following representation 

<^^/_i(z,t) = z U_____^ a(/) = 2. 
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We will end this section with results regarding the r-function representation of the second kind functions 
in the way we did in [12]. 



Lemma 7. The following identity 

n oo oo 

/\{Y1 = n A • • • A rn-i A ( XI '^n+i 

holds. 

Proof. Use induction in n. 

Theorem 11. Let ^ he a positive measure supported in T, then the following statements hold true 
1. The second kind functions have the following representation involving r-functions 



(110) 



□ 



r('+i)(t) 



R <\z 



C?{z,t) 



R- <\z 



Rl^ <\z 



a 



(0(t) 



-(0(t) 



< i?4 



< R. 



(Ill) 



<R-_\ 



2. For R- < \z\ < the Fourier series of the measure can be expressed in terms of r-functions in the 
following way 



2^r(0(t- [z-i]i) 



27rr(0(t+ [z]2) 
2^r(0(t+ [z]2) 



, a{l) = 1, 



(113) 



27rr(0(t- [z-i]i) 



, a{l) = 2, 



Proof. We will prove (jll2p only, and the proof of (jllip that goes analogously is left to the reader. The 
expression from Proposition H] can be arranged using the truncated columns of the moment matrix, that is 
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[I] .. 



/it X^'kx^''*)^d;u(t). Using this notation 

oo 



oo 



cJlAcflA-.-Acfl^A 5] z-'^^^^^cfs, 



3=1+1 

We define Z+^j := {j € Z+,a(j) = i}, i = 1,2 and notice that Z+ = Z+^i U The restrictions 

2^+12+15^^-12+2 °f mappings : Z+ i-+ N are bijections; hence, they have a well defined inverse, 

{v^)~^ and {v^)^^ . Therefore, 

oo 

r«(t)Ca(., t) = of A cfl A ... A cfl, A ^-cW^^..^^^,^^,^., 

3=0 



z-'^+('+^)c[;iAcflA...Acfl,A5]z-^cW 



where the last step requires the use of Lemma [7] with Cq , , . . . ;c|l,,c||^. Proceeding in a very similar 
way, we have 

oo 

dei{g^Ht))C^l{z,t) = r(^Ht)CP,{z,t) = c^ A A • • • A cfl, AY^cfx? 

3=1 

oo 

- Cq Ac, a • • • a c^_, a 2^ ^-"^ Cj- Oa(j),2 

oo 



A M A ... A A .^^(' + 2)-!'^ .iJ'] 

>_)-l(!^-(/+2+i)) 

j=0 



c[;UcfU...Acfl,Az-('+^)-i 



oo 



.-('+^)-ic[;iAcflA...Acfl,A^z^cf] 



(v_)-^{v.{l+2+3)) 



j=0 



with CqI , cf^l . . . , cf^j^l as adequate entries for Lemma [71 Finally, (I113p are obtained combining equations 
(jllip and (jll2p with Proposition [5l Theorem [9] and Theorem 1101 □ 

A comment on the convergence of the expressions in Theorem [TT] and it's proof is needed at this point. 
The main tool used in the proof is the series expansion of (1 — z)~^, that is convergent only for |z[ < 1 but 
can be analytically extended outside T. For instance, the r-function expression for \ is only strictly valid 
outside T. Nevertheless, it can be analytically extended inside the circle up to R- that is where the series for 
C2,i is convergent. As this extension is unique, we can talk about the analytically extended "Miwa-shifted" 
r-function. The same can be said about the other equalities involving r-functions; they are formally correct 
and convergent outside or inside T, but there is an analytical continuation for the shifted r-functions that 
converges where the Cauchy transforms do. The differences between expressions for the Cauchy transforms 
(jllip and (jll2p and their equivalent in the real line (e.g. [9j,|12j) is due to the existence of a positive and 
a negative part in the Laurent expansion around z = 0. The series expansion of (1 — z)~^ generates the 
positive part and the expansion of (1 — z^^)^^ gives a negative power series that generates the singular part 
of the Laurent expansion. 
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4.6 Bilinear equations 

For the derivation of a bilinear identity we proceed similarly as we did in [12] proving several lemmas. For 
the first one, let Wi,W2 be the wave matrices associated with the moment matrix g; so that, Wig = W2. 
Then, we have 

Lemma 8. The wave matrices associated with different times satisfy 

Wiit)Wi{tT^ = W2it)W2it'r\ (114) 

Proof. We consider simultaneously the following equations 

Wi{t)g = W2{t), 
Wi{t')g = W2{t'), 

and we get 

Wi{t)-^W2{t) = Wi{t')-^W2{t') = g, 
and the result becomes evident. □ 
Lemma 9. 1. For the vectors X^X* the following formulae hold 

Res,=o {X{X*V) = Res,=o = I, 

2. fT^ For any couple of semi-infinite matrices U and V we have 

UV = Res,=o{{Ux){V^X*V) (115) 
= Res,=o{{Uxl{V^xV) (116) 

We have the following 
Theorem 12. For any t,t' 

1. Wave functions satisfy 

Res,=o (^'i(z,t)(^'I(z,t'))^) =Res,=o [^2{z,t){'^*2{z,t')f 

2. OLPUC fulfill 

Res.=o (v.f^(.,t)c^«(.-\t')^-^i^.(^)e^-^(*^^^^'*^^^^^)) 

= Res.=oo {^['^ (z, t)^f {z-\t')z''F,{z) ^S^^^^-^^.^-') ) , (117) 

Proof. 1. First we notice that (jllSp and (jll6p can be written as 

UV = Res,=o {{Ux{z)){Vh*{z)Y 
= Res,=o (iUx*{z)){Vh{z))^ 



If we set in (fTTSll U = Wi{t) and V = Wi{t')-^ and in (fTTBIl we put U = W2{t) and V = W2{t')-^ 
attending to (|114p . recalling that ^1 = Wix, ^2 = W2X* and observing that = {W^^)^x* and 
^'2 = {W2^)^X we get the stated bilinear equation for the wave functions. 
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70 



2. We can substitute the expressions (l97|l and (198)1 to prove the second part of the result. 
We can reformulate this result using the residue theorem □ 

Corollary 4. If R- = and = oo and we take 70 and joo small zero-index cycles around z = and 
z = 00, respectively; then 

^^r\^,t)m)^^^Hz,t')dz= i ^^^\z,t){^l)^^\z,t')dz, (118) 
70 "^70 

^''^\z,t)^f{z~\t')e^T^^'^'^^''-*'2,^'')z~'F,{z)dz 

^'-^\z,t)^f{z-\t')e^T^^^'^^''~*'2,^-')z-^F,{z)dz, 

7cx) 

Alternatively the bilinear equation can he expressed using r-functions 

-'70 ^ ^ z 

= / r^^\,{t + [z-%)(T!iX'^{t - [z-%) + z'^+('+^)+-('+^)-V(;:2'^(t + [z],)) eEr^i(*^^-*2.>^- ^, 

if a(l)=l and 

i rfl,{t + [z],) ('+^)+V(5})(t + [z-\) + rfX'^it - [z],)) ^.^^^-^.^^^ ^ 

"'70 ^ ' z 

= I 4%,{t - [z]0(.-^-('-)-(^-)+V^S)(t - [z-%) + rtJ{t + Wi)) eS.-^fe.-y^^- 
J 70 / z 

if a(l)=2. 
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Appendices 

A Proofs 

Proof of Proposition [1] In this case g is positive definite and Hermitian, if we write 5*2 = hS2, where 

/ 1 (52)01 (52)02 . 

h = diag(/i07 ^1, • • • ) is a diagonal matrix and 82= 00 '^^'^^^'^ "' | , the uniqueness of the factorization 
implies that ^2 = {S^ ^)^ and /i^ E M and we get the stated result. 
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Proof of Proposition [3] Expressions (I16p . (I17p . (I19p and (j20p are obtained expressing the factorization 
problem as a system of equations. From (jl7p we deduce 

so that 

A;=0 ^ fc=0 

as stated in (fT8]l . To prove (f2T|) we consider 
SO that 

fc=0 ^^^^ yk=0 

that leads to (l2TD . 

Proof of Proposition [4] Using the definition for Ci, we have that 

k>l j=0 k>l j=0 



SO 



I 



j=0 



For the other set of functions we have 



I I 

k>l j=0 k>l j=0 

Comparing the expressions with those in Proposition[3]we see that they are formally identical, so we conclude 
the stated result. 



Proof of Proposition [B] Using the Fourier coefficients of the measure and the definition for the r„ ■ we 
obtain 



.(0 

LIVJ. tllC VJ-dllllUlVJil L\^L tllC J- 

k>0 k>0 k>0 



A;>0^ fc>0' 
fe>0 

k>-J{j) k<-J{j) 
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k>0 k>0 k>0 



k>0 

fc<J(i)+i fc>J(j)+i 



From here we obtain the rest of the expressions 



k>-J{j)+l k<~J{j)-l 
k<J{j)+l-l k>J{j)+l+l 

Proof of Proposition [8] From the definitions we have 



(Ci,i)t(z)vPi,i(z'] 


^=x\{z)Sl^SiXi{z') 


= {Xl)\z)xi{z') 


= Yz~'^~\zr 

n=0 


(C2,l)^(^)V'2,l(/) = 


xl{z)Sl{Sl)-\r{z') 


= {X\)\z)xi{z') 


oo 

= Y,z"^"\zr 

n=0 


{Cl,2)Hz)^,,,{z'] 


1 = X2{z)S^^Sixl{z') 


= {X2nz)x2{z') 


oo 
n=0 




X2{z)Sl{Slr\l{z') 


= {xl)Hz)xi{z') 


oo 
n=0 



which, after the study the region of convergence of the series involved, leads to the sated result. Then other 
identities derived from {x\)^ {z)x2{z') = {X2)'^ {z)xi{z') = 0. 

Proof of Proposition 1121 It can be deduced as follows 

J2k,2k+2 = {SiE2k,2k+2Si '^)2k,2k+2 = {Sl)2k,2k{Si ^)2k+2,2k+2 = Ij 

J2k,2k+1 = {SiE2k,2k+2Si ^)2k,2k+l = {Sl)2k,2k{S^ ^)2k+2,2k+l = — (5'l )2fc+2,2fe+l = —'^2k+2^ 

J2k,2k = {S2E2k+l,2k-lS2 ^)2k,2k = {S2)2k,2k+l{S2 ^)2k~l,2k = 

= —{S2)2k,2k{S2 ^)2k,2k+l{S2)2k+l,2k+l{S2 '^)2k-l,2k = ~'^^2k ^^2k+l^ 

J2k,2k-1 = (5'2-E'2fc+l,2fc-l5'2 "'^)2fc,2fc-l = (5'2)2fc,2fc+l(5'2 "^)2fc-l,2fc-l = 

= —{S2)2k,2k{S2 ^)2k,2k+l{S2)2k+l,2k+l{S2 ^)2k-l,2k-l = ~ P2k^^2k+1- 
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J2k+l,2k-l — (5'2-E2fc+l,2fc-l'S'2 ^)2k+l,2k-l — (5'2)2A:+1,2A:+1 (5'2 ^ ) 2k, 2k {S 2) 2k, 2k ^)2k~l,2k-l 

_ 2 2 

— P2k+lP2k^ 

J2k+l,2k = (5*2-^2/^+1, 2fc-l5'2 ^)2k+l,2k = (5'2)2fc+l,2fe+l (5*2 "^)2fc-l,2fc 

= (5'2)2fc+l,2fc+l(5'2 ^)2k,2k{S2)2k,2k{S2 ^)2k-l,2k 

_ 2 -(2) 

— P2k+l'^2k ' 

J2k+l,2k+l = (5'l-E2fc,2fc+2'S'i )2fe+l,2fe+l = {Sl)2k+l,2k{Si "'")2fc+2,2fc+l = " ('S'l )2fc+l,2fe (^l )2fc+2,2fc+l 



-(2) (1) 
"'^2fc+l'^2fc+2' 



<^2A:+l,2A;+2 — ('S'l-E2fc,2fc+2'S'i )2fc+l,2fc+2 — {Sl)2k+l,2kiSi ^)2k+2,2k+2 — '^2k+l- 

Proof of Proposition 1141 For k = the result comes from the definition of p^. For k = 1,2 we have to 
use the truncated recursion relations, 

_1 (0) (1) -(2) (0) (0) (2) (1) (1) (1) (0) 

^-1^(1) _ ^(3) _ ^(2)^(2) _ (1) (2) (1) _ 2-(2) (0) (1) _ -(2) (2) _ (1) (2) (1) 2 (0) 

multiplying by z and integrating we obtain ho = hi — a^^^ J ipf'^ zdfi, and J (^['^^zd/i = — a^^^/io, from 

where we have ho = hi + hQ. Now multiplying by z^^ and integrating we obtain = a['^^/i2 — 

02^'' a^^'* / ip^-^^ z~^dfi + Pif ipi'^z^^dp, f ip^i'^z'^dfi = —p\a2^hQ and / iJi^ z'^dp = —Oi^hQ leading to 

= /i2 + a2^a'2^hi — hi. 

The other k > 2 can be proved by induction. For odd case we multiply (j34p by z^^^ to obtain 

n (1)^ (l)-(2) f (2fc) fc+l- , 2 (1) /" (2fe-l) fc+ln ,2 2 /" (2fc-2) fc+ln 

= "2fc^2fc+l-a2fc "2fc+l / V^i ^ d// + P2fca^fc_l / f/?! ^ d/Z + P2fc-lP2fc / V^i ^ d/x 

Jt Jt Jt 

then multiplying by z'^ the recurrence expressions (|3ip and (j32p for ZLp'^^\ zip^^ -zv^^^*^ and integrating 
we substitute, to get 

(\ U I (1) -(2) u (1) -(2) I 2 7 

= /l2fc+l + a2fc+ia2fc+i/l2fc - a2fc-ia2fc-l'i2fc - P2fc-1^2fc, 

from where using the induction principle the result is proven. Should we want to obtain the rest of the 
equations (those with even k\ then it is necessary to multiply by z~^~^ the odd recurrence relations for 
z^p^^^^"^^ (p2|) and use the same procedure. 

Proof from Proposition 1171 From the definition we have 

1-1 /-I 



\ ifW(n,z')i^W(^,n)dM(n) := / Y.^f {z')pf {u)Y^pf {u)^f {z)dp{u) 

T J'^3=^ k=0 

^rf)(z')|:#(^-) / ^?{u)^^^\u)dp{u) 



j=0 k=0 

l-l l-l 



^^?iz')J2^^^\-z)6,,, = kV\z,z'). 

j=0 k=0 
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Alternate proof for Theorem [3] Due to Gerardo Ariznabarreta. 

From the block factorization problem g^^^ = {S^^)~^ ■ ^2 is upper-triangular, we can write + 

then using the definition for Lp'^^ and using the previous formula, we get 

3=0 r,j=0 k>0 

r,j=0 

In addition, we can express the formula for ip^^~^^^ in the following way as 

j=0 
j=0 

r,j=0 fc>0 

+ (x('+') - (5m,o ••• 5m,/-i) (5t'^)-'x"l), 

from where we obtain (if / is odd) ipi\i{z) = ipf~^^\z) — ai-i-iipi\z) . 

Proof of Lemma [1] If we denote by IfW = Yl\=o (^^^ projection over the first / components) we find 
i^W(z,/) = (nH$2(^))^nW$i(z') = <l>2{z)^U^^^^i{z') = x{z)^ s^^uV^Sixiz'). 

From the block factorization g^^^ = {sf^)^^S2^ and its inverse (g^^^)^^ = {S2^)^^sf^ we can get an expression 
for using only finite size matrices, that is 

K^Hz,z') = xizM^S2'U^'^S^U^\iz') 

= x^H^)\s]^Y'sf\v}(^z') 

Proof of Lemma [2] The symmetry of g in pop can be expressed using the block structure 
T -- 
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g[l] 


g[l,>l] 







With this block structure we get 

fV]g[i] + '):[i,>i]g[>i,i] = gWxW +^['.>']t[^'''1, 

or equivalently, recalling the Gaussian factorization, we arrive to 

(^W)-ItH - rV^gV])-'^ = (^W)-i(c,['.>'lTfe'''] - T^^'>^]g[>^fi)(gV])-'^. 
We have also the equations 

that leads to 

{z' - -z-^)K^%z,z') = X^Hz)\g^^^)-^z'x^%z') - -z-\^\z)Hg^^Y\^\z') 
= xW(^)n(5W)-^TW-TH(5M)-^)xW(^')+ 

Proof of Proposition [19] If a(/) = 1 then z''-'^'-^ipf\z) is a monic polynomial of degree + ~ 

1, while when a{l) = 2 then is a monic polynomial of degree + i^+(0 ~ 1- The 

orthogonality relations for z'^-^''^ipf\z) and z'^+^''^'^(pf\z^^) can be obtained from ()65p 



T 



z'^-«4'\(z)z-'=dM^) = 0, A: = 0,...,|z7(0!-1, a(0 = 1, 



(119) 



that means that 



z-W^«(z) = P|,^(,)l_i(z), a{l) = l, 

z'^^('^-'^g{z~')=P\,^l^\_,(z), a{l) = 2, 

recalling that |z7(/)| — 1 = / we get the desired result. 



Proof of Proposition j 

1. The shift operators defined fulfill 

^n,lXn{z) = zUfi^iXniz), ^n,2Xn{z) = z''^Ufi^2Xniz) , 

^H,lXn{z) = (Z~^n^,i - ^0,0A"^)Xn(^), ^]i,2Xn{z) = (2^2 " (A^ )"+ )x^(z) , 

that means that 

(AS,1 + AT + En+,n+{A^T+)Xn{z) = ZXn{z), (AT, + A.H,2 + ^0,oA"+ )Xn(z) = z'^Xniz), 
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from there it follows that 



5n(Ali + A.j,2 + ^o,oA"+)^ 



<7n(A^,l+Ajo + ^n+,nH-(A'^)"+; 



2. With the definitions 



the use of Proposition [23] leads easily to 



J, 



n,2 



n,l 



(120) 



The matrix x has n_ + 1 diagonals over the main diagonal and the matrix 2 has n+ + 1 diagonals 
under the main diagonal (in both computations we have excluded the main diagonal itself), so both 
JniiJn2 have 72+ + n_ + 3 diagonal band. 



Proof of Proposition 1281 We have calculated previously J = Li and using the same method L2 can be 
calculated, both are five-diagonal matrices given by 

\ 



L,=J 







1 

















pI 




-(2) 

a\ 




















-p14^ 


(1) -(2) 
-"3 «2 




1 














pIpI 


2 -(2) 
^302 




-(2) 

Q3 




















2 (1) 


{l)-(2) 


-4" 


1 














pIpI 


2 -(2) 
P5«4 




-(2) 
"5 






V ■■ 



-(2) 


1 




















2 -(2) 

-pfa'2 


{l)-{2) 


-(2) 

—03 
















P1P2 


P2(^i^ 


(l)-{2) 

-"2 "3 






















2 -(2) 
-P3«4 


(l)-(2) 

— Q4 


-(2) 
















pIpI 


2 (1) 
Pi«3 


(l)-{2) 






















2 -(2) 




-(2) 

— Qy 


1 
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as we are looking only for time flows associated to tu and t2i then Bi i = (Li)+ and B2.1 = (-L2)-. Using 
([93|) in the matrix elements {Li)k.k+i for /c > we obtain 



da 



(1) 

2k+2 



l)2k,2k+l 



dt 



da 



11 
(1) 

2k+2 



dt. 



da. 



21 

(2) 

2fe+l 



dtu 

d{Li)2k,2k+l 
)2k+l,2k+2 



-[Bl,l,Li\2k,2k+l - a2fc+3i^ ~"2fc+2"2fc+2J 



-[S2,l,Li]2fc,2fc+l - "21+1(1 - a2fc+2"2fc+2 



„(1) ^,(2) 



dt 



11 



5t 



11 



[D r 1 -(2)n (1) -(2) 

[Bi^l,Li\2k+l,2k+2 = (1 ~ "2fc+l" 



2fc+l^ 



^•^L+l _ 5(Li)2A;+l,2A:+2 



dt 



21 



5t 



21 



[^2,1,^ 



lj2fc+l,2fc+2 



"2A:+2V-^ "2fc+l"2A;+lV 



and now looking at {Li)2k,2k and (^2)2^+1,2^+1 for A; > we obtain the rest of the equations 



^(«2fc«2fc+l) _ d{Li)2k,2k 



dt 



11 



9(«2fc«2fc+l) 



dt 



21 



^("2fc+l"2fc+2) 
5tll 

^('^2fc+l'^2/c+2) 



5tii 

d{Li)2k,2k _ _ 
5t21 

^(-^2)2fc+l,2fc+l 
C^(-^^2)2fc+l,2fc+l 



dt 



21 



21 



[-Bl,l, -i^l]2fc,2fc ^ 
[-82,1, il]2fc,2fe ^ 
= — [^1,1, -^^2]2A:+1,2A:+1 
= — [-62,1, i2]2A:+l,2A:+l 



da 



(1) 

2fc+l 



dt 



da. 



11 
(1) 

2fc+l 



(1) /I _ (1) -(2) N 
'^2fc+2V-'^ '^2A;+l'^2fc+lJ 



dt 



da. 



21 

(2) 

2A;+2 



4^(1 - a 



2fc+l"2fc+l> 



dt 



da. 



11 

(2) 

2fc+2 



dt 



21 



"2fe+iv-^ "2fc+2"2fc+2'' 



-(2) _ (1) -(2) N 

~^2k+^\^ '^2k+2^2k+2> 



considering all the equations we obtain ()99p . 

Proof of Proposition [T8] First let us look to (jSOp . Using Definition 1161 it can be expressed as 



i,j=0 



l-l 



^-i^(x('+")(z)det5W - Yl gi^U-^y'-'Mffx'^'H 



j=0 i=0 

that is the expansion of ([SU|) . Using the same idea with 

I , i 



we arrive at the expansion of (j53p taking the complex conjugate. Both (j5ip and ()52p can be proved using 
the same ideas. 
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